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The  problem  of  determining  robust,  low  order  controllers  which  achieve  stability  and 
disturbance  attenuation  is  considered.  This  is  an  important  problem  for  the  control  of  very 
high  order  structures  where  a  full  order  controller  is  inappropriate. 

To  address  this  problem,  an  approach  using  the  Frobenius-Hankel  (FH)  norm  is  devel¬ 
oped.  The  FH  norm  is  shown  to  lead  to  attractive  robustness  and  performance  properties. 
A  parameter  optimization  method  is  developed  which  provides  an  iterative  method  for  de¬ 
termining  the  FH  optimal  parameters  of  a  dynamic  system  using  a  gradient  approach. 

Two  distinct  controller  design  methods  are  presented  which  make  use  of  the  FH  optimiza¬ 
tion  procedure.  The  first  approach  determines  a  controller  which  minimizes  the  FH  norm  of 
the  closed  loop  system.  While  the  solution  method  is  iterative,  the  procedure  proves  to  be 
straightforward  to  apply. 

The  second  approach  uses  projective  controls  as  a  first  step  in  the  design.  This  allows 
the  designer  to  take  advantage  of  the  attractive  features  of  projective  controls.  Design 

parameters  in  the  controllers  are  then  determined  by  solving  the  FH  optimization  problem 

/  .  .  -  " 

called  the  ''‘auxiliary  minimization  problem.'” 

Examples  of  each  approach  are  given,  including  a  full  design  problem  for  the  control  of 
a  flexible  structure  using  the  projective  controls  approach.  The  resulting  improvement  in 
the  disturbance  attenuation  of  the  system  using  only  second  order  controllers  points  to  the 
effectiveness  of  this  design  procedure.  - 
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CHAPTER  1 
INTRODUCTION 

The  design  of  low  order  controllers  to  achieve  stability,  attenuation  of  disturbances 
and  robustness  to  both  modeling  errors  and  parameter  variation  is  an  important  design 
problem  in  the  control  of  many  complex  systems.  This  is  particularly  true  in  the  control 
of  large  flexible  structures  where  factors  related  to  cost,  implementation,  maintenance 
and  reliability  may  make  a  collection  of  low  order,  decentralized  controllers  the  preferred 
configuration. 

The  goal  is  thus  to  determine  a  robust,  low  order  controller  such  that  the  closed  loop 
system  is  stable  and  has  good  performance.  The  crucial  constraint  dealt  with  in  this  work 
is  on  the  order  of  the  controller  which,  while  not  necessarily  predetermined,  is  to  be  of  order 
significantly  lower  than  the  order  of  the  system.  However,  the  problem  of  determining  the 
best  such  controller  under  these  restrictions  is  an  open  problem  at  present. 

The  general  control  problem  considered  in  this  work  is  illustrated  in  Figure  1.1:  Given 
the  system  P(s),  determine  a  low  order,  possibly  decentralized  controller  C(s)  so  that 
the  closed  loop  system  has  good  performance  and  robustness  properties.  The  constraints 
on  the  order  of  the  controller  are  usually  dictated  by  limitations  on  the  computational 
power  of  the  processor  which  implements  the  controller.  Decentralized  controllers  often 
are  nictated  when  physical  obstacles  prohibit  the  availability  of  all  the  measurements  at 
one  central  location. 

Classical  control,  which  relies  on  root  locus  and  Bode  plots,  uses  a  one-loop-at-a-time 
method  to  design  controllers.  This  is  quite  inadequate  for  MIMO  systems  of  high  order. 
Thus,  for  complex  systems,  ;t  is  common  to  design  controllers  based  on  optimization 


Figure  1.1:  General  control  problem. 


9 


techniques.  This  approach  allows  complex  controllers  to  be  designed  using  a  single  criterion 
of  optimality.  The  difficulty,  however,  lies  in  choosing  the  correct  criterion  and  computing 
the  resulting  optimal  control. 

Considerable  success  has  been  realized  in  computing  optimal  controllers  for  norms  such 
as  H?  (LQG)  [12]  and  Hoo  [5].  In  general,  closed  form  solutions  for  the  optimal  controller 
can  be  found  in  terms  of  the  solution  of  Riccati  equations.  However,  the  optimal  controller 
in  general  is  a  centralized  controller  of  order  at  least  as  high  as  the  plant.  In  many  cases 
this  is  not  practical.  This  is  particularly  true  in  the  control  of  flexible  systems  where  the 
plant’s  inherent  high  order  makes  the  use  of  such  controllers  unacceptable. 

Unfortunately,  the  computation  of  low  order  and  decentralized  optimal  controllers  has 
not  been  as  successful  as  the  full  order  centralized  case.  The  computations  of  such  con¬ 
trollers  are  much  harder  to  solve  since  at  present,  there  are  no  direct  methods  of  computing 
the  optimal  controller  as  in  the  full  order  case  and  most  methods  of  computing  optimal 
controllers  are  either  iterative  or  yield  suboptimal  controllers. 

A  parametric  approach  to  the  problem  is  to  parameterize  the  controller  and  then  use 
standard  minimization  techniques  to  optimize  the  parameters  of  the  controller.  Such  meth¬ 
ods  have  been  extensively  used  in  finding  low  order  Hi  optimal  controllers  [13,31,11,14]. 
Difficulties  with  this  approach  include  local  minima  and  slow  convergence.  Also,  this 
method  is  difficult  to  apply  in  the  H^  case  due  to  the  difficulty  in  directly  computing  the 
h i^j  norm. 

A  method  for  computing  an  optimal  controller  of  fixed  order  which  does  not  involve  pa¬ 
rameterizing  the  controller  has  been  shown  for  the  Hi  optimal  case  [9]  and  for  the  bounded 
Hoo  norm  case  [3].  This  approach  has  the  advantage  that  the  solution  method  is  indepen¬ 
dent  of  the  controller  realization.  However,  the  computation  of  the  optimal  controller 
involves  the  solution  of  two  coupled  Riccati  equations  and  two  Lyapunov  equations.  At 
present,  an  efficient  method  of  solving  these  matrix  equations  is  unavailable  which  limits 
its  application  for  truly  high  order  systems. 

Because  of  the  difficulty  in  computing  optimal  controllers,  suboptimal  methods  are 
quite  common  due  to  their  relative  ease  of  computation.  However,  in  general,  these  meth¬ 
ods  will  not  always  produce  near  optimal  solutions  and  may  in  fact  produce  controllers 
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which  destabilize  the  closed  loop  system. 

One  common  suboptimal  method  involves  the  use  of  model  reduction.  Model  reduction 
can  be  used  to  compute  low  order  controllers  in  two  distinct  ways.  One  way  is  to  approxi¬ 
mate  the  plant  model  with  a  reduced  order  plant  model  and  then  to  use  a  full  order  design 
method  to  compute  the  resulting  lower  order  controller.  The  other  way  is  to  first  compute 
the  full  order  controller  and  then  approximate  the  full  order  controller  with  a  lower  order 
controller.  These  approaches  are  neither  optimal  nor  guaranteed  to  be  stable.  However, 
they  can  be  significantly  simpler  to  compute  than  the  optimal  solution.  For  an  overview 
of  these  methods,  see  [1]. 

As  can  be  seen,  no  one  method  exists  which  produces  optimal  or  near  optimal  low 
order  controllers  in  an  efficient  computational  manner.  In  particular,  no  method  exists 
to  compute  low  order  Hoo  optimal  controllers  which  are  attractive  due  to  the  robustness 
properties  of  the  norm. 

This  work  develops  the  use  of  the  Frobenius-Hankel  norm  as  an  optimization  criterion 
for  the  design  of  low  order  and  decentralized  controllers.  This  norm  is  established  as  a  useful 
norm  for  this  task  through  the  examination  of  its  analytical  and  numerical  properties.  The 
FH  optimization  approach  is  used  in  two  different  ways:  First,  to  compute  FH  optimal 
controllers  directly,  and  second,  by  using  it  to  solve  the  subproblem  of  determining  the 
free  parameters  of  projective  controllers.  In  both  cases,  the  centralized  and  decentralized 
controller  cases  are  each  examined. 

The  organization  of  this  work  is  as  follows:  First,  a  review  of  system  theory  and  norms 
is  given.  Then  the  Frobenius-Hankel  norm  is  introduced  and  its  properties  are  explored.  In 
particular,  its  relationship  to  the  7^^  and  the  7i2  norm  establish  its  robustness  properties. 
The  following  chapter  develops  necessary  conditions  for  an  optimal  FH  norm  solution  and 
develops  numerical  methods  to  solve  the  optimization  problem. 

The  next  chapter  investigates  the  computation  of  low  order  optimal  controllers.  Pre¬ 
viously,  no  techniques  existed  for  the  computation  of  low  order,  'H<X)  optimal  controllers. 
In  this  chapter,  necessary  conditions  for  the  solution  of  the  optimal  controller  are 
derived.  However,  computational  problems  exist  in  solving  these  conditions.  This  leads 
to  the  consideration  of  the  use  of  FH  optimization  as  a  means  of  determining  low  order 
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controllers.  The  results  of  the  previous  chapter  are  applied  to  the  optimal  control  problem 
to  determine  necessary  conditions  for  an  optimal  solution.  The  optimization  procedure 
of  the  previous  chapter  is  demonstrated  by  calculating  an  FH  optimal  controller  for  an 
example  problem. 

The  use  of  projective  controls  is  then  investigated  as  a  method  of  computing  low  order 
controllers.  Although  this  is  not  a  new  technique,  there  axe  a  number  of  differences  between 
the  published  work  and  the  presentation  in  this  dissertation.  The  perspective  here,  is  an 
input/output  approach  rather  than  the  traditional  state  space  approach.  In  particular,  a 
number  of  input /output  relationships  are  developed  which  give  insight  into  the  problem. 
A  simpler  formulation  of  the  proper  controller  case  is  given  which  eliminates  the  need  to 
transform  the  plant  into  a  specific  form.  The  FH  optimization  method  is  then  applied  to 
the  problem  of  determining  the  free  parameters  of  the  projective  controllers.  The  following 
chapter  examines  the  case  of  decentralized  projective  controllers. 

In  order  to  demonstrate  this  approach,  a  design  example  using  projective  controls  and 
FH  minimization  to  determine  the  free  parameters  is  presented.  The  control  problem 
co  isists  of  controlling  the  flexible  modes  of  a  cruciform  satellite  structure.  Using  a  40- 
state  analysis  model,  two  decentralized,  second  order  controllers  are  determined  which 
attenuate  the  effect  of  the  disturbance. 


1.1  Notation 


ft",  Cn 
at,  a* 

Tr  (A) 

A;(A) 

<7  ,(A),a(A),£(A) 

n2 


I|G(s)||2 

||G(*)||oo 

I|G(-)IIf 


n-dimensional  real  and  complex  Euclidean  spaces 
transpose  and  complex  conjugate  transpose  of  A  6  TZnxm 
trace  of  A  €  7£n*n  ,  Def.  7 

ith,  largest  and  smallest  eigenvalue  of  A  €  7ZnXn  resp. 
ith,  largest  and  smallest  singular  value  of  A  6  7ZnXm 
space  of  all  stable,  strictly  proper  transfer  functions 
space  of  all  stable,  proper  transfer  functions 
H2  norm  of  G(s)  G  7^2,  Def.  4 
Hoo  norm  of  G(s)  €  Hoc ,  Def.  5 
Frobenius-Hankel  norm  of  G(s)  €  H2,  Def.  6 


ith  Hankel  singular  value  of  G(s ),  Def.  3 
Dirac  delta  function  (impulse) 


o 


S(t) 


A 

B 

C 

D 

Cl’}, 


transfer  function,  Def.  1 

Laplace  and  Fourier  transform  operators,  resp. 
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CHAPTER  2 
BACKGROUND 

2.1  Representation  of  Dynamic  Systems 


A  linear,  time-invariant  (LTI)  system  is  usually  represented  in  one  of  two  equivalent 
ways.  The  first  involves  expressing  the  input-output  relationship  of  the  system  in  terms 
of  first  order  differential  equations.  This  is  often  called  the  state  space  representation  of 
the  system  because  of  the  use  of  the  state  variable.  The  other  way  to  represent  the  system 
is  through  the  Laplace  transform  of  the  impulse  response  of  the  system.  This  is  called  the 
transfer  function  of  the  system. 

While  both  are  equivalent  representations  of  an  LTI  system,  each  has  its  own  strengths 
and  weaknesses  as  a  system  representation.  In  particular,  the  transfer  function  repre¬ 
sentation  works  well  to  represent  the  output  of  the  system  as  a  function  of  the  system 
inputs.  Thus,  it  emphasizes  the  input-output  nature  of  the  system.  However,  the  transfer 
function  of  the  system  is  defined  in  terms  of  polynomials  of  a  complex  variable.  Thus,  it 
poses  computational  problems  for  use  in  design  and  analysis.  The  state  space  approach 
is  useful  because  it  can  easily  be  defined  in  terms  of  real  matrices.  This  makes  the  state 
space  system  useful  for  computational  purposes.  Thus,  it  follows  that  it  would  be  useful 
to  be  able  to  use  the  compact  notation  of  the  transfer  function  while  retaining  the  state 
space  representation  for  computational  uses.  The  following  notation  is  introduced  for  this 
purpose. 


Definition  1  Given  the  state  space  representation  of  the  LTI  system 


i(t)  =  Ax(t)  +  Bw(t )  (2.1) 

z(t)  =  Cx(t)  +  Dw(t )  (2.2) 


where  x  €  7ln  is  the  state,  w  €  is  the  input,  z  6  TV  is  the  output,  the  transfer  function 
of  the  system  is  denoted 


A 

B 

C 

D 

=  C(sI-  A)~'B  +  D. 


(2-3) 


Thus,  through  the  use  of  the  expression  (2.3),  the  transfer  function  of  the  system  is 
represented  in  terms  of  a  state  space  representation  using  a  compact  notation. 


2.1.1  Series  connection 


Gi(s) 


G2(s) 


Figure  2.1:  Series  connection. 

Given  two  systems  Gx(s)  and  G2(s) 

a  Bx  a  A2  B2 

Gx(s)  ±  _J - L  ,  g2(s)  ^  — - — - 


c>  Dr 


C2  D2 


connected  together  in  series  as  shown  in  Figure  2.1,  the  resulting  system  can  be  represented 
as 

-  * 

A2  B2Ci  B2D\ 

G2(s)Gx{s)=  0  Ax  Bx  ■  (2.5) 

C2  D2Cx  D2Dx 

Note  that  this  representation  is  not  necessarily  minimal. 

2.1.2  Parallel  connection 


Gx(s) 


G2(s) 


Figure  2.2:  Parallel  connection. 

Given  two  systems  Gi(s)  and  G2(s)  connected  in  parallel  as  in  Figure  2.2,  the  resulting 
system  can  be  represented  as 


4i 

0 

Bx 

0 

A2 

b2 

c, 

c2 

D\  4-  D2 
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Note  that  this  representation  is  not  necessarily  minimal. 

2.1.3  State  space  transformation 

A  state  space  transformation  has  no  effect  on  the  transfer  function  of  a  system.  Thus, 
given  the  state  space  transformation  x  =  Tx  where  T  is  a  nonsingular  matrix, 


GM  = 


A 

B 

'  T-1AT 

T~'B 

C 

D 

CT 

D 

(2.7) 


This  points  out  an  important  fact  about  state  space  representations  of  a  system.  That 
is,  they  are  not  unique  representations  of  a  transfer  function. 

2.1.4  Minimal  representation 

The  minimal  representation  of  a  system  is  one  such  that  it  has  no  uncontrollable  or 
unobservable  states.  If  a  system  contains  uncontrollable  and  unobservable  states,  it  can 
be  reduced  to  its  minimal  representation  by  removing  the  uncontrollable  and  unobservable 
states. 

This  can  be  easily  done  if  the  system  has  the  following  canonical  form: 


<?oo  = 


Ac  a 

Al2 

Al3 

An 

Bco 

0 

Aco 

A23 

A24 

Bco 

0 

0 

A(0 

Am 

0 

0 

0 

0 

Aco 

0 

0 

C CO 

Cco 

0 

D 

(2.8) 


This  system  is  equivalent  to  the  following  system 

G(s)  = 


■ 

- 

Aco 

Bco 

_  Cco 

D 

(2.9) 


Note  that  AcS  represents  dynamics  which  are  controllable  but  not  observable,  Aco  represents 
dynamics  which  are  both  controllable  and  observable,  ASo  represents  dynamics  which  are 
observable  but  not  controllable,  and  Aes  represents  dynamics  which  axe  neither  controllable 
nor  observable. 


9 


P(s) 


C(s) 


Figure  2.3:  Feedback  connection. 

2.2  Feedback  Systems 


Given  the  feedback  system  system  shown  in  Figure  2.3  where  the  plant  Pis )  has  the 


state  space  representation 

x  =  Ax 
z  =  Hx 
y  =  Cx 


■f  Bu  +  Gw 

(2.10) 

+  Eu 

(2.11) 

+  Dw 

(2.12) 

and  x  G  TV  is  the  state,  w  €  TV  is  the  disturbance  input,  u  €  TV1  is  the  controlled  input, 
z  €  TV  is  the  controlled  output,  and  y  €  IV  is  the  measured  output,  the  transfer  function 
P(s)  such  that 


z(s) 

=  P(s) 

w(s) 

.  y . 

.  u ^  . 

(2.13) 


is  given  by 


A 

G  B 

P(s)  = 

H 

0  E 

c 

D  0 

(2.14) 


Note  that  the  control  system  described  here  assumes  that  the  controlled  outputs  do  not 
contain  a  direct  disturbance  term  and  that  the  measured  outputs  do  not  contain  a  direct 
control  term.  These  assumptions,  however,  are  not  overly  restrictive  since  they  represent 
common  assumptions  for  control  problems.  Also,  the  omission  of  these  terms  does  not 
form  a  fundamental  restriction,  but  rather  one  made  for  notational  simplification. 

If  the  system  P(s )  is  controlled  by 


w>  = 


Ac  Bc 

Cc  Dc 


(2.15) 
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2.3  Hankel  Singular  Values 

The  Hankel  singular  values  and  norms  defined  on  Hankel  singular  values  of  a  linear 
system  have  been  shown  to  have  useful  applications  in  model  reduction  [23,7],  in  providing 
bounds  for  the  TC0 0  norm  [7],  and  in  the  design  of  low  order  controllers  [24,19,27]. 

Definition  2  The  controllability  and  observability  grammians  of  a  system  G{s)  6 
7f2  where 

G(s)  = 

are  defined  as,  reap. 

P=  r  eMBBTeATtdt  (2.19) 

Jo 

Q  =  H  eATtCTCeAtdt.  (2.20) 

Jo 

The  grammians  of  the  system  are  the  unique  positive  definite  solutions  to  the  Lyapunov 
equations 

AP  +  PAT  +  BBT  =  0  (2.21) 

AtQ  +  QA  +  CTC  =  0.  (2.22) 

These  equations  form  an  efficient  approach  to  solving  for  the  grammians. 

Definition  3  The  Hankel  singular  values  of  G(s)  6  7t2  are  given  by 

»i(G(*))  =  IAi(/’Q)]'/’  (2.23) 

where  P  and  Q  are  the  controllability  and  observability  grammians  of  G(s). 
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2.4  Norms  of  Dynamic  Systems 

The  two  most  common  norms  used  for  dynamic  systems  are  the  W2  norm  and  the 
norm.  These  norms  are  reviewed  here  for  reference  in  later  chapters. 

Definition  4  The  H2  norm  of  G{s)  €  H2  w  defined  as 

hgwiii  =  [t  {i  jT  cow)"Go«)<fc>}]  .  (2.24) 

The  W2  norm  is  often  called  the  “quadratic  norm”  due  to  its  interpretation  as  the  integral 
of  the  square  of  the  impulse  response. 

Theorem  1  Given  G(s)  €  H2  and  the  impulse  response  g(t)  =  [£?(«)] 

l|G(s)||l  =  Tt  r  (2.25) 

JO 

The  'H.2  norm  can  be  computed  from  a  state  space  representation  of  the  system. 
Theorem  2  Given  G(s )  6  TL2  and  the  controllability  and  observability  grammians  P  and 

Q, 

||G(s)||l  =  Tr  {PCTC}  =  Tr  {QBBT}.  (2.26) 

The  'H.oo  norm  is  often  considered  a  “worst-case”  norm  and  as  such  has  applications  to 
robustness  issues  since  it  can  be  used  to  bound  the  performance  of  a  system. 

Definition  5  The  Tioo  norm  of  G(s)  €  'H00  is  defined  as 

IIGMIU  =  8up*[C(,'«)].  (2.27) 

LJ 

The  Hoc  norm  is  the  equivalent  to  the  maximum  gain  of  the  system  as  shown  in  the 
following  theorem. 

Theorem  3  Given  G(s )  €  'Hoo 

IIGMIU  =  sup  f±  (2.28) 

u/*0  ||W>||2 

where  z(t)  is  the  response  of  G(s)  to  the  input  w(t). 
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This  theorem  is  an  important  robustness  result  since  it  states  that  if  a  bound  is  known  on 
the  norm  of  the  input,  a  bound  can  be  calculated  for  the  norm  of  the  output. 

The  computation  of  the  Hoo  norm  is  a  difficult  problem  since  it  cannot  be  computed 
directly.  However,  an  upper  bound  for  the  H<x,  norm  can  be  computed  as  follows  [4,5]: 


Theorem  4  Given  G(s)  £  Hoo  such  that 

G{s)  = 


A 

B 

C 

D 

(2.29) 


and  7  >  <x(D).  Then  HG^)]^  <  7  iff  H  has  no  imaginary  eigenvalues  where 


H  = 


A  +  BR'1DTC  y~2BR~lBT 
-CTS~XC  -AT  -CTDR~1BT 


(2.30) 


and  R*  I  -  7  ~2DTD  and  S  =  I  -  7  ~2DDT. 

For  the  case  where  the  system  is  strictly  proper,  the  theorem  can  be  simplified  as  follows: 
Theorem  5  Given  G(s )  €  H?  such  that 


<?«  = 


A 

B 

C 

0 

(2.31) 


Then  IjG^^loo  <  7  iff  H  has  no  imaginary  eigenvalues  where 


H  = 


A  7  ~2BBT 
-CTC  -AT 


(2.32) 


Equivalently,  ||G(s)||oo  <7  if  and  only  if  there  exists  an  X  >  0  which  satisfies  the  Riccati 
equation 

AtX  +  XA  +  7  ~2XBBtX  +  CTC  =  0.  (2.33) 


Thus,  the  Hoo  norm  of  a  system  can  be  computed  in  an  iterative  way  by  applying 
a  search  algorithm  to  the  problem  of  finding  the  smallest  7  which  satisfies  the  bound 
->jjG(s)||oo  <  7 •  This  7  is  the  Hoo  norm  of  the  system. 
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CHAPTER  3 

THE  FROBENIUS-HANKEL  NORM 


Recently,  Medanic  and  Perkins  introduced  the  Frobenius-Hankel  (FH)  norm  for  the 
design  of  control  systems  [17].  The  FH  norm  is  defined  as  the  Frobenius  norm  on  the 
Hankel  singular  values.  The  motivation  for  its  use  stems  from  its  relationship  to  more 
widely  known  norms  such  as  Hi  and  'Hoa  and  its  good  computational  properties  which 
make  it  suitable  for  use  in  optimization  procedures. 

In  this  chapter,  the  Frobenius-Hankel  norm  is  defined  and  its  properties  explored.  In 
particular,  both  a  time-domain  and  frequency- domain  physical  interpretation  will  made 
of  the  meaning  of  the  FH  norm  and  a  simple  computational  method  will  be  shown  for 
calculating  the  FH  norm.  The  FH  norm  will  also  be  directly  related  to  both  the  7f2  and 
‘Htx  norms. 

In  the  following  chapters  the  FH  norm  will  be  used  as  the  basis  for  a  parameter  op¬ 
timization  problem  and  applied  to  a  model  reduction  problem  and  the  low  order  optimal 
controller  problem. 


Theorem  6  Given  the  system  G(s )  £  K2  and  the  controllability  and  observability  gram- 
mians  P  and  Q  resp.,  then 

||G(s)||2F=Tr{PQ}.  (3.2) 


Proof:  From  Definition  6, 


||G(s)||2F  =  TrE2 


(3.3) 
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where  £  =  dia.g(a\...o n).  Since  there  exists  T  nonsingular  such  that  T  T PT  1  =  £  and 
TQTt  =  £  [7] 

||G(s)|j£  =  Tr  (T'TTT)£(rT_1)£  (3.4) 

=  Tr  (Tt£T)(T-1  ET~t)  (3.5) 

=  Tr  PQ.  (3.6) 


Note  that  this  expression  involves  the  solution  of  two  Lyapunov  equations  but  avoids 
computation  of  the  Hankel  singular  values.  This  expression  for  the  FH  norm  is  important 
because  the  trace  function  is  more  well-behaved  than  the  singular  value  function. 


3.1.1  Time-domain  property 'S 

A  time-domain  interpretation  of  the  Frobenius-Hankel  norm  is  as  follows: 

Theorem  7  Given  the  system  <7(s)  €  W2  and  the  impulse  response  of  the  system  g(t), 
then 

IIGMIIJ.  =  TV  rtg(t)Tg(t)dt.  (3.7) 

JO 

Proof:  From  Theorem  6, 

IIGMIIJ.  =  Tr  PQ.  (3,8) 


From  Definition  3, 


=  lim  Tr 

[T  eAtBBTeATt  dt 

f  eATrCTCeAr  dr 

T  — *oo 

Jo 

Jo 

which  is  equivalent  to 

Tr  PQ  =  Hm  Tr  £  £  [CeA(t+T)B\  [CeA(t+r)B\T  dt  dr. 
Let  g(r)  =  CeMB , 

rT  rT 

Tr  PQ  =  lim  Tr  /  /  g(t  +  r)T g(t  +  r)  dt  dr 

T—oo  Jo  Jo 

Tr  PQ  =  lim  Tr  f  f  g(t)T g(t)  dt  dr 

T—oo  Jo  JT 


(3.9) 


(3.10) 


(3.11) 

(312) 
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Let  H(r)  =  /tT+t  g(t)T g(t)  dt. 

TtPQ  =  lim  Tr  F  H(t)  dr. 

T —co  Jo 

Integrating  by  parts, 

Tr  PQ  =  lim  Tr  |tf(r)r|J  -  F  t  dH(r) 

T—oo  [  Jo 

Tr  1  Q  =  lim  Tr  f  (T  -  t)g(t  +  T)T g(t  +  T)  +  tg(t)T  g{t)  dt 

T-+00  JO 

In  the  limit  as  T  — *  oo,  g(t  +  T)  — >  0,  thus 

Tr  PQ  =  Tr  f°  tg(t)T g(t)  dt. 

Jo 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


Note  the  similarity  of  this  expression  to  the  one  given  in  Theorem  1.  In  fact,  a  comparison 
of  the  two  expressions  shows  that  the  FH  norm  is  in  fact  equivalent  to  a  time- weighted  H.2 
norm.  The  weighting  factor  t  acts  to  weight  relatively  more  heavily  the  impulse  response 
at  later  times  than  at  earlier  times.  Thus,  large  initial  responses  are  weighted  less,  but 
lightly  damped  dynamics  are  weighted  more. 


3.1.2  Frequency-domain  properties 

Theorem  8  Given  the  system  G(s )  €  Ti-i  and  the  frequency  response  of  the  system  G(juj)  = 


G(5)l,=n,>  ihen 


\\G(s)\\2f  =  ^Tr  ^pG(j^r  du. 

Proof:  Applying  Parseval’s  Theorem  to  (3.7)  yields 

iicMiiJ-  =  ^Tr 


=  r  j 

2tt  J- oo  j  \  du>  t 


G(ju>Y  du> 


(3.17) 


(3.18) 

(3.19) 


Note  that  the  this  expression  shows  that  there  are  two  components  to  the  FH  norm.  The 
first  component  comes  from  the  term  in  the  above  expression  which  implies  that  the 

FH  norm  weights  the  “flatness”  of  the  transfer  functions.  The  second  component  is  due 
to  the  G(jw)  term  in  the  above  expression  which  implies  that  the  FH  norm  also  weights 
the  magnitude  of  the  transfer  function. 
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3.2  Relationships  with  Other  Norms 

Although  a  physical  interpretation  has  been  given  for  the  Frobenius-Hankel  norm, 
relating  it  to  familiar  norms  is  useful  in  providing  further  motivation. 

The  Frobenius-Hankel  norm  can  be  related  to  the  Hoo  norm  through  the  Hankel  singular 
values  of  the  system.  In  the  following  two  theorems,  expressions  are  shown  which  relate 
the  two  norms  using  the  Hankel  singular  values. 

This  first  result  defines  an  interval  in  which  both  the  Hoo  and  FH  norm  must  lie.  Thus, 
results  which  apply  to  the  Hoo  norm  can  be  applied  to  the  FH  norm  through  this  result. 

Theorem  9  Given  G(s)  E  Hi, 


n 


»(GM)<||G(*)|U<22>(C(a)) 

1=1 

(3.20) 

»(G(a))  <  ||G(a)||r  <  y;«r,(G(j)). 

1=1 

(3.21) 

Proof:  The  first  expression  has  been  shown  in  [6]  and  [7].  To  show  the  second  expression, 

begin  first  with  the  upper  bound. 

CM  *  (M  UN 

n  n 

(3.22) 

= 

i=i  j=i 

n  n  n 

(3.23) 

=  + 

fc=l  •»!  j= 1 

'*) 

n 

(3.24) 

>  ■£  4  =  IIGWIIJ.. 

k=  1 

The  lower  bound  follows  directly  from 

(3.25) 

=  E^,2  >  o. 

1=1 

(3.26) 

■ 

Using  the  results  of  the  previous  theorem,  the  H0 o  norm  can 

be  bounded  in  terms  of 

the  FH  norm.  This  gives  a  measure  of  the  “closeness”  of  the  Hoo 

norm  to  the  FH  norm. 

Theorem  10  Given  G(s )  E  Hi, 

4=I|GWIIf  <  IIGWIW  <  2>/S||GW||r. 

y/  n 


(3.27) 
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Proof:  To  show  the  lower  bound,  first  note  that 

J2  <  rid2. 
i=i 

Then  it  follows  from  (3.20)  that 

IIGMIU  >  s  >  -4||CWIIf 

y/n 

which  proves  the  lower  bound.  To  show  the  upper  bound,  note  that 

f>,  =£T£ 

i=l 

1  1  ...  1  j  and  ET  = 


where  ET  = 


(T\  <72 


(3.28) 


(3.29) 


(3.30) 

The  Frobenius  norm  of  a  matrix 


is  defined  as 


||X||Fr  =  [TrATA],/2. 

Since  the  Frobenius  norm  is  a  valid  matrix  norm 

(3.31) 

I|J3TS||f,  <  ||ET||Fr||E||F. 

(3.32) 

=  [Tr£T£]1/2[TrETE]1/2 

(3.33) 

=  V"!!G(s)||f. 

(3.34) 

However,  since  ETE  is  a  positive  valued  scalar, 

\\ETE\\ft  =  ETZ.  (3.35) 

Thus, 

±*><VZ\\G(s)\\f.  (3.36) 

i=i 

Substituting  into  (3.20)  yields  the  upper  bound.  ■ 

These  expressions  establish  the  robustness  properties  of  the  FH  norm  through  the 
norm.  Thus,  since  the  computation  of  the  FH  norm  is  much  simpler,  one  rationale  for  the 
use  of  the  FH  norm  is  established. 

The  Frobenius-Hankel  norm  can  also  be  related  to  the  sensitivity  of  the  norm  to  a 
shift  of  the  eigenvalues  of  the  system  along  the  reed  axis.  Such  a  shift  in  the  eigenvalues 
can  be  represented  by  letting  A  be  given  by 

A(a)  =  A0  +  aI.  (3.37) 
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Thus,  A(A)  =  A(A0)  +  q  and  the  sensitivity  of  the  Hi  norm  to  a  shift  in  the  eigenvalues 
along  the  real  axis  is  given  by  ^||G(s)||2.  The  following  theorem  relates  this  sensitivity  to 
the  FH  norm. 


Theorem  11  Let  the  eigenvalues  of  the  system  G(s )  €  hii  be  given  by  \t  =  A,  4-  a,  then 


Proof:  Let 


where 


Then 


where  PQ  =  &  satisfies 


Let  Q  satisfy 


i-iiGWi  -2KGWIU. 

A  '  Or=0 

(3.38) 

J  t  ||G(s)||*  =  Tr  PCTC 

(3.39) 

AP  +  PAt  +  BBt  =  0. 

(3.40) 

^  =  Tr  PaCTC 
da 

(3.41) 

APa  +  PaAT  +  2P  =  0. 

(3.42) 

QA  +  AtQ  +  CTC  =  0. 

(3.43) 

Then,  using  the  properties  of  the  trace,  it  can  be  shown  that 


Thus 


Tr  PaCTC  =  2Tr  PQ. 


^  =  2Tr  PQ. 
da 


(3.44) 


(3.45) 


This  expression  establishes  the  FH  norm  as  a  sensitivity  measure  for  the  H?  norm  with 
respect  to  a  change  in  the  relative  stability  of  the  system.  This  is  an  important  relationship 
since  in  flexible  structures  problems,  often  the  damping  of  the  system  is  not  known  well. 
Also,  a  change  in  the  relative  stability  of  a  lightly  damped  structure  is  the  variation  most 
likely  to  cause  instability. 
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CHAPTER  4 

FROBENIUS-HANKEL  NORM  OPTIMIZATION 


In  this  chapter  a  common  framework  for  solving  optimal  FH  norm  problems  will  be 
presented.  In  later  chapters,  specific  problems  will  be  solved  under  this  framework. 

The  general  problem  to  be  solved  is:  given  a  parameterized  system,  find  the  set  of 
parameters  which  minimize  the  FH  norm  of  the  system.  The  first  section  involves  finding 
necessary  conditions  for  such  a  set  of  optimal  parameters.  Solutions  to  these  necessary 
conditions  are  explored  in  the  following  sections.  The  final  section  applies  the  general  FH 
optimization  problem  to  the  problem  of  FH  optimal  model  reduction. 

4.1  The  FH  Optimization  Problem 

The  FH  optimization  problem  is:  Let  the  parameterized  system  G(s)  £  7i.2  be  given  by 


where  £  is  the  parameter  matrix.  The  FH  norm  of  G(s)  can  be  computed  as 

J  =  ||<?(«)l£  =  Tv{PQ}  (4.2) 

where  P  and  Q  satisfy 

U  =  AP  +  PAt  +  BBt  (4.3) 

0  =  AtQ  +  QA  +  CTC.  (4.4) 

The  optimization  problem  is  thus,  find  £  such  that  the  criterion  (4.2)  is  minimized  subject 
to  the  constraints  (4.3)-(4.4). 

4.2  General  Solution  Methods 

This  constrained  optimization  problem  can  be  converted  to  an  unconstrained  optimiza¬ 
tion  problem  using  Lagrange  multipliers.  The  augmented  criterion  is  given  by 

J  =  Tr  {PQ  +  L(AP  +  PAt  +  BBT)  +  M{ATQ  +  QA  +  CTC)|. 


(4.5) 
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Using  this  approach,  necessary  conditions  for  an  optimal  solution  are 
dJ 


=  AtL  +  LA  +  Q  =  0 
=  AM  +  MAt  +  P  =  0 


dP 
dJ_ 

dQ 

-  =  AP  +  PAt  +  BBt  =  0 
oL 

^L=AtQ  +  QA  +  CtC  =  0 
oM 

^  =  ^Tr  {2 At(LP  +  QM)  +  BTLB  +  CTCM}  =  0. 


(4.6) 

(4.7) 

(4.8) 

(4.9) 
(4.10) 


In  general  these  equations  can  not  be  solved  for  the  optimal  if  directly.  However,  iterative 
methods  may  be  applied  to  this  problem. 


4.2.1  A  steepest  descent  approach 


A  steepest  descent  approach  to  this  problem  is  to  find  the  direction  of  steepest  descent 
and  to  take  a  step  in  that  direction.  The  direction  of  steepest  descent  is  in  the  direction 
of  the  gradient  wrt  £.  The  gradient  of  J  wrt  to  £  is  given  by 


dJ  dJ  dJ  dP  dJ  dQ  dJ  dL  dJ  dM 
d£  ~  d£  +  dP  d£  +  dQ  d£  +  dL  d£  +  dM  d£  ‘ 

If  P,Q,L,M  satisfy  (4.6)-(4.9),  then 

d£  ”  5£- 


The  parameter  update  is  given  by 


C.-+1  =  Zi  - 


(4.11) 


(4.12) 


(4.13) 


This  method  was  is  illustrated  in  Figure  4.1  and  used  in  [27]  to  solve  the  FH  optimization 
subproblem. 


4.2.2  Riccati  equation  approach 

The  Riccati  approach  [24]  uses  Riccati  equations  instead  of  Lyapunov  equations.  The 
Riccati  equations  are  constructed  so  that  the  iterative  solution  converges  to  the  solution 
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1.  Select  f0  so  that  A(£0)  is  stable. 

2.  Let  i  =  1. 

3.  Solve  Eqs.  (4.6)-(4.9l  for  L.  M,  P  and  Q. 

4.  Calculate  £,+1  from  Eqs.  (4.12)  and  (4.13). 

5.  If  the  parameters  have  not  converged,  let  i  =  i  +  1  and  go  to  2. 


Figure  4.1:  Steepest  descent  algorithm. 


of  the  Lyapunov  equations.  The  advantage  of  this  approach  is  that  the  equations  can 
be  solved  for  all  systems,  stable  or  unstable.  Thus,  an  initial  stabilizing  controller  is  not 
needed.  However,  the  drawback  to  this  approach  lies  in  the  fast  that  Riccati  equations  are 
computationally  much  more  time-consuming  to  solve  than  Lyapunov  equations. 


The  iterative  equations  are 

.A7'.Lt-+i  +  Ll+\A  —  Li+iRL,+i  4*  L{RL ,  4-  Q  =  0  (4.14) 

AMi+1  4-  Mt+1AT  -  Mi+lRMi+x  +  M,RM,  4-  P  =  0  (4.15) 

APt+l  +  Pi+lAT  -  Pi+1RPi+l  4-  PRP,  +  BBt  =  0  (4.16) 

ATQi+ 1  4-  Q,+i-A  —  Qi+iRQi+\  4-  Q,RQ,  4-  CTC  =  0  (4.17) 


where  is  the  solution  of 


(4. IS) 


Note  that  if  this  iterative  algorithm  converges,  it  converges  to  the  solution  of  the  corre¬ 
sponding  Lyapunov  equations.  Thus,  a  solution  of  these  equations  satisfies  the  necessary 
conditions  for  an  optimal  solution. 


4.3  FH  Optimal  Model  Reduction 


The  model  reduction  problem  is  an  important  one  since  it  is  often  necessary  to  approxi¬ 
mate  a  high  order  model  with  a  lower  order  model.  The  use  of  a  lower  order  approximation 
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Figure  4.2:  Model  reduction  problem. 

makes  design  and  analysis  simpler  although  with  a  less  accurate  model. 
The  model  reduction  problem  is  as  follows:  Given  an  n-th  order  system 


G(s)  = 


find  a  &-th  order  approximation 


G(s)  = 


A  B 
C  D 


A  B 
C  D 


that  minimizes  the  FH  norm  of  the  error  ||G(s)  —  G(.s)||f\ 
The  error  system  as  illustrated  in  Figure  4.2  is 


G(s)  -  G{s)  = 


Ae  Be 
Ce  De 


A  0 


=  0  A 


-C  D-D 


(4.19) 


(4.20) 


(4.21) 


Note  that  the  parameters  of  the  error  system  are  simply  the  matrices  that  define  the 


reduced  order  model.  Thus, 


A  B 
C  D 


(4.22) 


In  order  to  apply  the  results  for  FH  optimization,  the  system  to  be  optimized  must  be 
strictly  proper.  Thus,  for  the  case  of  model  reduction,  the  error  system  must  be  strictly 
proper.  This  is  satisfied  if  and  only  if  De  =  0.  Thus, 


D  =  D. 


(4.23) 


In  order  to  simplify  the  solution  of  this  problem,  note  that  the  matrices  which  define 
the  error  system  can  be  written  as  a  linear  function  of  the  parameters  £.  Thus. 


Ae  = 


0  A 


Ik  0  A  B 


0  0  C  D 


Ik  0 


=  -40  +  TiT 


(4.24) 
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h 

• 

0 

A 

B 

0 

Be  = 

+ 

c 

B 

L  J 

0 

0 

D 

.  Im  . 

r 

r  i 

A  B 

i - 

o 

0  c 

+ 

o 

1 

• 

L 

C  D 

1 

o 

o 

De  =  D  + 


r 

A  B 

o 

1 

b" 

L 

C  D 

—  B0  +  T£B  i 

=  C0  +  C^T 

=  D0  +  Cy^Bi- 


(4.25) 

(4.26) 

(4.27) 


This  greatly  simplifies  the  necessary  conditions  since  the  parameters  axe  given  by  a  single 
matrix. 

The  necessary  conditions  from  (4.6)-(4.9)  are 


| Z  =  A.P  +  PAT,  +  B.BTt  =  0 

(4.28) 

§jfi  =  A'Q  +  QA-  +  C?C'  =  0 

(4.29) 

^  =  a*l  +  la.  +  q  =  0 

(4.30) 

=  A  'M  +  MAT,  +  P  =  0. 

(4.31) 

Since  the  parameters  £  are  given  by  (4.22),  (4.10)  becomes 

C\  j 

—  =  2  [T(LP  +  QM)T  +  TLBeBj  +  C\C,MT  =  0.  (4.32) 


The  FH  optimal  model  reduction  could  thus  be  solved  by  applying  one  of  the  methods 
discussed  in  the  previous  sections. 
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CHAPTER  5 

THE  H0 o  AND  FH  OPTIMAL  CONTROLLER 

In  this  chapter,  two  methods  of  computing  low  order  optimal  controllers  are  developed. 
First,  necessary  conditions  for  an  Hoo  optimal  controller  are  derived.  These  represent  new 
results  in  the  determination  of  low  order  H oo  optimal  controllers.  However,  computational 
problems  make  this  approach  difficult. 

The  second  method  developed  is  FH  optimal  control.  This  method  is  chosen  due  to  its 
robustness  properties  as  shown  in  Chapter  3  and  its  good  numerical  properties. 

5.1  The  Hoo  Optimal  Output  Feedback  Controller 

The  ur  j  of  the  Hoo  norm  is  motivated  by  robustness  concerns.  The  H0 0  norm  of  a 
system  describes  the  maximum  gain  from  input  to  output  of  the  system.  Thus,  it  gives 
the  “worst-case”  performance  of  the  system.  This  is  very  useful  particularly  when  very 
little  is  know  about  the  input  signal.  The  Hoo  norm  also  is  useful  since  it  can  be  used  to 
guarantee  stability  of  the  system  despite  uncertainty  in  the  model  of  the  system. 

This  problem  has  been  solved  for  the  state  feedback  and  full  order  feedback  cases  [5]. 
Also,  sufficient  conditions  for  computing  a  low  order  controller  which  satisfies  an  upper 
bound  on  the  Hoo  norm  of  the  system  have  been  shown  [3,2].  The  goal  of  this  chapter  is 
to  determine  necessary  conditions  for  the  solution  of  the  Hoo  optimal  control  problem. 

The  Hoo  optimal  control  problem  is  as  follows:  Given  the  plant  P{s),  determine  the 
controller  C(s)  which  minimizes  the  Hoo  norm  of  the  closed  loop  system  G{s). 


Let  the  plant  have  the  form 


The  controller  is  assumed  to  be  a  static  controller.  Thus, 

C(s)  =  K.  (5.3) 

The  closed  loop  system  then  is 


5.1.1  State  feedback  controllers 

The  problem  of  determining  the  Hen  optimal  controller  in  the  case  of  state  feedback 
(C  =  I)  has  been  shown  in  [5].  The  results  are  repeated  here  for  comparison  with  the 
following  section. 

Theorem  12  The  Hen  bound  ||G(s)||oo  <71/  and  only  if  there  exists  M  >  0  such  that 

AtM  +  MA  +  Q  +  7~2MGGtM  -  MBR~lBTM  =  0.  (5.5) 

If  such  an  M  exists,  then  one  controller  which  satisfies  the  Hoo  bound  is 

C{s)  =  K  (5.6) 

where 

K  =  -R~iBTM.  (5.7) 

Note  that  the  conditions  in  the  theorem  are  necessary  and  sufficient.  Thus,  the  controller 
which  corresponds  to  the  smallest  7  which  satisfies  the  conditions  of  the  theorem  is  the  Hen 
optimal  controller.  Thus,  the  state  feedback  Hen  optimal  control  problem  can  be  solved 
by  using  a  search  algorithm  on  7  to  determine  the  minimum. 

5.1.2  Output  feedback  controllers 

The  case  for  the  case  of  output  feedback  controllers  is  not  as  simple.  This  thesis 
represents  the  first  work  on  Hen  optimal  output  feedback  controllers.  Previous  work  focused 
only  on  sufficient  conditions  for  an  Hen  bound  to  hold. 

From  Theorem  5,  the  bound  ||G(s)||oo  <  7  holds  iff  there  exists  M  >  0  such  that 


AtM  +  MA  +  7~2MBBtM  +  CTC  =  0. 


(5.8) 


Thus,  the  1-Loo  optimal  control  problem  reduces  to  the  constrained  minimization  prob- 


inf  J  =  72  (J 

subject  to  (5.8)  and  M  >  0. 

Using  Lagrange  multipliers,  the  equivalent  unconstrained  minimization  problem  is 


inf  J  =  72  +  Tr  L(AT M  +  MA  +  i~2MBBtM  +  CTC)  +  £  ntXx(M). 

K  « 


(5.10) 


The  necessary  conditions  for  an  optimal  solution  axe  the  Kuhn*Tucker  conditions  given 


riiXi(M)  =0,  n,  >  0 


(5.11) 


i-n  =  (i  +  -,-2BBtU)L  +  LI  A  +  +  XNX~'  =  0  (5.12) 

oM 


=  AtM  +  MA  +  7 -2MBBtM  +  CTC  =  0 
:  =  2RBtMLCt  +  2KCLC7  =  0 


~^  =  1~  7~4Tr  LMGGtM  =  0 
a- y 

where  N  =  diag{rii}  and  X  satisfies  MX  =  XA.  with  A  =  diag{\i(M )}. 


(5.13) 

(5.14) 

(5.15) 


The  controller  parameter  K  can  be  eliminated  by  introducing  the  projection  matrix  r. 
From  (5.14), 

K  =  -R-'BTMLCT(CLCT)~l.  (5.16) 

Then  (5.13)  becomes 

AtM  +  MA  +  7 ~2MGGtM  +  Q-  MBR-'BtM  +  tImBR~1BtMtx  =  0  (5.17) 


where 


r  =  LCt(CLCt)-'C 

Ti  =  I  —  T. 


(5.18) 

(5.19) 


From  (5.15),  an  expression  for  7  may  be  derived: 


72  =  [Tr  LMGGtM}$. 


(5.20) 
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However,  this  expression  is  not  useful  for  determining  the  optimal  7  since  L  cannot  be 
determined  uniquely.  Note  that  if  L  =  L0  is  a  solution  to  (5.12),  then  L  =  aL0  is  also 
a  solution.  Thus,  the  optimal  controller  is  the  controller  associated  with  the  minimum  7 
such  that  the  necessary  conditions  have  a  solution. 

Note  that  from  this  expression  one  can  conclude  that  if  the  minimum  7  is  non-zero, 
then  L  /  0.  But  from  (5.12),  this  is  only  possible  if  A  +  ~i~2MGGT M  has  purely  imaginary 
eigenvalues  or  M  becomes  indefinite.  Thus,  the  minimal  7  can  be  found  by  reducing  the 
value  of  7  until  one  of  the  stopping  conditions  is  violated.  The  stopping  conditions  are 

1.  A  +  7 ~2MGGT M  has  no  purely  imaginary  eigenvalues 

2.  M  >  0 

In  the  state  feedback  case,  C  is  full  rank  and  thus  r  =  I.  This  decouples  (5.13)  and 
(5.12).  Thus  (5.13)  becomes 

AtM  +  MA  +  ~f-2MGGTM  +  Q-  MBR‘1BTM  =  0.  (5.21) 

Note  that  this  agrees  with  the  results  for  the  state  feedback  case  given  in  [5]. 

One  possible  algorithm  for  solving  this  problem  is  given  in  Figure  5.1.  This  algorithm 
has  proved  to  have  slow  convergence  properties  which  become  more  acute  as  the  minimum 
is  approached.  However,  efficient  algorithms  for  the  solution  of  these  equations  will  require 
additional  investigation  beyond  the  scope  of  this  work. 

5.2  The  FH  Optimal  Controller 

This  section  derives  the  necessary  conditions  for  an  FH  optimal  controller.  Note  that 
the  assumption  of  the  use  of  a  time-invariant  controller  is  in  itself  a  restriction.  Appendix  B 
shows  that  the  unconstrained  controller  is  a  time- varying  controller.  Thus,  the  assumption 
of  a  time-invariant  controller  is  in  itself  a  constraint  and  results  in  a  higher  FH  norm  than 
in  the  unconstrained  case. 

The  FH  optimal  control  problem  is  as  follows:  Given  the  plant  P(& )  determine  the 
controller  C(s )  which  minimizes  the  norm  of  the  closed  loop  system  G(s). 
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1.  Let  t  =  /,  7  =  oo. 

2.  Find  smallest  7  such  that  (5.17)  has  a  solution. 

3.  Solve  for  L  from  (5.12). 

4.  Iteratively  solve  (5.17)  and  (5.18)-(5.19)  for  r  and  M. 

5.  Go  to  2. 

Figure  5.1:  An  Hoo  optimal  controller  algorithm 


Let  the  plant  have  the  form 


A 

G 

B 

H 

0 

E 

C 

D 

0 

The  controller  is  assumed  to  be  a  p-th  order  dynamic  controller.  Thus, 

^  ,  ’  Ac  Bc  ' 

C(s)= - . 

Cc  Dc 

The  closed-loop  system  is 


G(a)  = 


A  B 
C  D 


A  +  BDcC  BCc  G  +  BDCD 
BcC  Ac  BCD 


(5.22) 


(5.23) 


(5.24) 


L  1  J  [  H  +  EDcC  ECc  I  EDCD  J 
The  parameters  of  the  closed  loop  system  are  the  matrices  that  define  the  controller.  Thus, 


_  Ac  Bc 

’  Cc  Dc  ' 

In  order  to  simplify  the  notation,  the  following  expressions  are  used: 


(5.25) 


A  =  A  +  BiC 


(5.26) 
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B  =  G  +  BiD 
C  =  H  +  E£C 
D  =  E£D 


(5.27) 

(5.28) 

(5.29) 


where 


0 

0 

.  r  iP  o 

Ip  o' 

A  = 

0 

A 

B=\  P 
[  0  B 

C  = 

0  c 

0 

r  i 

0 

D  = 

D 

o 

II 

Kl 

G  = 

G 

H  = 

[  0  H 

. 

(5.30) 


In  order  for  the  closed-loop  system  to  be  strictly  proper  it  is  required  that  D  =  0. 
Thus,  one  condition  for  the  optimal  solution  is 


ED,D  =  0. 


The  necessary  conditions  for  an  optimal  solution  from  (4.6)-(4.9)  are 

=  AP  +  PAt  +  BBt  =  0 
oL 


dJ 

dM 


=  AtQ  +  QA  +  ClC  =  0 


dJ  ~t 

—  =  atl  +  la  +  q  =  o 

J 

-—  =  AM  +  MAt  +P  =  0 
oQ 


Equation  (4.10)  becomes 

dj 


—  =  2  [Bt{LP  +  QM)CT  +  BtL(G  +  B£D)Dt  +  ET(H  +  E^C)MCT}. 


(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 

(5.36) 


If,  in  addition,  the  controller  is  non-dynamic,  i.e.,  C(s)  =  Dc,  then  the  closed  loop 
system  is 


A  +  BDcC 

G  ' 

A 

B 

H  +  EDcC 

0 

C 

D 

G(s)  = 


The  necessary  conditions  for  an  optimal  control  reduce  to  (5.32)-(5.35)  an<^ 


(5.37) 


dJ 

dDc 


=  2  [Bt(LP  +  QM)Ct  +  BtMG  +  BDcD)Dt  +  Et{H  +  EDcC)MCt}.  (5.38) 


Satisfying  the  condition  EDCD  =  0  leads  to  three  special  cases: 
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Figure  5.2:  Plant  with  decentralized  controllers 

1.  Noise-Free  Measurements  (D  =  0) 

2.  Cheap  Control  ( E  =  0) 

3.  Strictly  Proper  Controller  ( Dc  —  0) 

If  the  measured  outputs  of  the  plant  do  not  contain  a  direct  feedthrough  term  for  the 
disturbance,  then  D  =  0.  This  is  referred  to  as  the  “noise-free”  case  since  the  measured 
outputs  are  considered  to  be  uncorrupted  by  any  disturbance.  If  the  controlled  outputs 
of  the  plant  do  not  include  a  control  term,  then  E  =  0.  This  is  referred  to  as  the  “cheap 
control”  case  since  the  application  of  large  amounts  of  control  energy  is  not  reflected  in 
the  controlled  outputs.  If  the  measured  outputs  include  a  non-singular  noise  term  ( DDT ) 
and  the  controlled  outputs  include  a  non-singular  control  term  (ET E),  then  Dc  =  0.  In 
this  case,  the  resulting  controller  is  strictly  proper. 

Of  course  (5.31)  can  be  satisfied  in  a  combination  of  these  three  cases.  This  would 
imply  that  some  channels  of  the  control  would  be  strictly  proper,  some  noise-free  and 
some  with  cheap  control. 

5.3  FH  Optimal  Decentralized  Controllers 

The  decentralized  control  problem  shown  in  Figure  5.2  arises  when  all  the  controls  and 
measurements  are  not  available  at  one  central  location.  Thus,  each  separate  channel  has 
a  separate  controller  which  works  with  the  available  measurements  and  controls. 

The  FH  optimal  decentralized  control  problem  is:  Given  a  system  P(s)  whose  measure¬ 
ments  and  controls  have  been  decentralized  into  /-channels  controlled  by  /  decentralized, 
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dynamic  controllers,  find  the  l  controllers  C,(s)  which  minimize  the  Frobenius-Hankel  norm 
of  the  closed  loop  system  G(s). 

Assume  P(s)  has  the  form 


P(s)  = 


A 

G 

Bt  . 

..  Bt 

H 

0 

Ex  . 

■  E, 

c , 

0  . 

..  0 

C, 

Di 

0  . 

..  0 

(5.39) 


and  the  /  controllers  are  given  by 


C,(s)  = 


’  Ad 

Bd  ' 

ca 

Dci 

(5.40) 


The  solution  to  this  problem  is  to  reformulate  it  in  terms  of  the  centralized  control  prob¬ 
lem.  Thus,  define  the  centralized  controller  C(s)  in  terms  of  the  decentralized  controllers 

Ci(s). 

'  Ci(s)  0 


C(s)  = 


0 


Ci(s) 


'  Ac 

Bc 

Cc 

L>c 

(5.41) 


where 


Ac  =  diag(Acx, . 

Ac/ ) 

(5.42) 

Bc  =  diag(Bcx, 

... ,Bcl ) 

(5.43) 

Gc  ~  diag(Cc\,. 

..,Cc) 

(5.44) 

Dc  =  diag(Dcl , 

(5.45) 

Also,  let 


Cx 

Dx 

C  = 

,  D  = 

* 

c, 

r>i 

B  =  B\  ...  B , 


(5.46) 


(5.47) 
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E=  \  Ei  ...  Ei].  (5.48) 

. 

Thu?  the  decentralized  optimal  control  problem  is  identical  to  the  centralized  case 
except  that  the  parameters  of  the  controller  axe  constrained.  Necessary  conditions  for  an 
optimal  solution  are  thus  (5.32)-(5.35)  with  the  additional  parameter  constraints  given  by 

(5.42M5.45). 

5.4  Example 

In  order  to  demonstrate  the  computation  of  an  FH  optimal  controller,  an  example 
control  problem  is  presented. 

5.4.1  Problem  definition 

The  optimization  problem  is:  Given  the  plant 


P(s)  = 


'  -0.4335 

-0.0118 

-0.9231 

-0.4643 

0.8854 

-0.7382  ' 

-0.9160 

-0.5185 

-0.4110 

-0.0779 

0.1747 

1.5473 

-0.0414 

-0.6085 

-0.7507 

-0.8901 

-1.4939 

0.8204 

-0.4828 

-0.0916 

-0.2014 

-0.9215 

-1.1423 

-1.5361 

0.9782 

1.9938 

-0.8140 

-0.8819 

0 

-1.5443 

0.1821 

0.3387 

1.6250 

1.0326 

0 

0 

(5.49) 


determine  a  second  order,  proper,  stabilizing  controller  C(s)  which  minimizes  the  FH  norm 
of  the  closed  loop  system. 

5.4.2  Controller  optimization 

To  start  the  algorithm,  an  initial,  non-zero,  stabilizing  controller  is  needed.  Such  a 
controller  is  given  by 


-1  0 
C(s)  =0-1 
0.0846  -0.1728 


0  -0.0118 
-1  -0.0555 
728  0 


(5.50) 
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Figure  5.3:  Iteration  history  of  FH  norm  of  the  system. 


Using  the  Matlab  routine  given  in  Appendix  C,  the  optimal  controller  was  computed. 
Figure  5.3  shows  the  FH  norm  at  each  iteration  in  the  algorithm  along  with  the  Hoc  bounds 
of  the  system  given  in  (3.20).  In  this  example,  the  FH  norm  was  reduced  from  its  ini  tied 
value  of  84.9  down  to  5.6.  and  the  upper  bound  on  the  Hoo  norm  reduced  from  240  down 


to  20.  The  optimal  controller  is 

determined  to  be 

-0.7966  -0.2337 

-0.4563 

C(s)  = 

-0.2502  -0.7149 

0.5331 

0.4515  -0.5542 

0.1822 

The  frequency  response  of  both  the  open  loop  and  closed  loop  systems  is  given  in  Figure 
5.4.  As  can  be  seen,  the  Hoo  norm  has  been  reduced  from  55  to  5.  Note  that  while  the 
Hoo  bounds  in  this  case  were  a  bit  conservative,  FH  optimization  was  quite  effective  in 
reducing  the  Hoo  norm  of  the  system. 
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CHAPTER  6 

PROJECTIVE  CONTROLLERS 

The  projective  controls  method  [15,16]  is  a  method  for  designing  low  order  controllers 
for  higher  order  systems.  This  approach  has  been  studied  for  a  variety  of  approaches 
[8,21,20,22,30,29].  The  essence  of  the  method  is  that  it  restricts  the  class  of  low  order 
controllers  considered  to  the  subclass  of  controllers  that  retain  the  dominant  eigenvalues 
and  eigenstructure  of  a  reference  state  feedback  controlled  system.  Thus,  given  a  reference 
system  with  eigenvalues  given  by 

Are/  =  ArJ  U  Ar2  (6-1) 

and  the  associated  eigenstructure  span{Xri}  and  span{Xri},  the  closed  loop  system  under 
projective  controls  has  eigenvalues 

A proj  —  Ari  U  AP2  (6.2) 

and  retains  the  eigensubspace  span{Xri},  where  Arl  contains  the  retained  eigenvalues  from 
(6.1)  and  Ap2  contains  the  remaining,  i.e.,  residual  eigenvalues  of  the  closed  loop  system. 
The  family  of  controllers  which  achieve  this  axe  conveniently  parameterized  by  a  single  free 
parameter  matrix  of  dimension  p  x  r  where  r  is  the  number  of  available  measurements  and 
p  is  the  order  of  the  controller. 

The  retention  of  the  eigenstructure  as  well  as  the  eigenvalues  of  the  reference  system  is 
particularily  important  in  disturbance  attenuation  since  the  orientation  of  the  eigenspaces 
influences  the  transfer  function  and  system  zeros.  Thus,  if  the  state  feedback  reference 
solution  is  determined  to  provide  disturbance  attenuation,  this  is  reflected  in  the  eigen¬ 
structure,  and  the  retention  of  the  eigenstructure  together  with  the  eigenvalues  is  an  added 
incentive  to  the  use  of  projective  controls  in  disturbance  attenuation. 

This  effect  of  projective  controls  can  be  characterized  through  the  transfer  function  of 
the  reference  system.  Let  the  transfer  function  from  the  disturbance  input  to  the  regulated 
ouput  for  the  state  feedback  controller  reference  system  be  given  by  Gr(s).  This  transfer 
function  can  be  reduced  to  the  series  connection 


Gr(s)  =  Grl(s)Gri(s) 


(6.3) 
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where  G>i(s)  contains  the  dynamics  of  the  reference  system  to  be  retained  and  Gr 2(s) 
contains  the  remaining  dynamics.  The  corresponding  transfer  function  under  projective 
controls  can  be  reduced  to  the  series  connection 

Gp(s)  =  Gn(s)Gp2(s)  (6.4) 

where  Gri(s )  is  called  the  retained  subsystem  since  it  contains  the  dynamics  retained 
from  the  reference  system  and  Gp2(s)  is  called  the  residual  subsystem.  The  order  of 
the  retained  subsystem  is  determined  by  the  class  of  controller  chosen.  Once  the  order 
and  dynamics  to  be  retained  sire  chosen,  the  computation  of  the  projective  controller  is 
straightforward.  Three  classes  of  controllers  axe  considered  here:  static,  proper  and  strictly 
proper  controllers. 

The  reference  system  is  determined  by  a  state  feedback  controller  which  is  chosen  for 
its  desirable  properties.  Many  algorithms  exist  for  designing  state  feedback  controllers, 
thus  the  projective  controls  method  is  suitable  for  use  in  combination  with  many  types 
of  synthesis  methods.  Moreover,  once  the  state  feedback  controller  is  determined,  the 
projective  controller  is  easily  computed.  However,  stability  and  performance  of  the  residual 
dynamics  are  not  guaranteed. 

In  the  case  of  dynamic  controllers,  free  parameters  exist  which  may  be  used  to  im¬ 
prove  the  performance  of  the  residual  subsystem.  Section  6.5  addresses  the  problem  of 
determining  the  free  parameters. 

In  this  chapter,  expressions  axe  shown  for  the  projective  controllers,  the  closed  loop 
eigenvalues  and  the  error  between  the  reference  system  and  the  projective  system. 

6.1  The  Reference  System 

A  state  feedback  controller  u  =  K0x  applied  to  the  plant  yields  the  reference  system 


where  F  =  A  +  BK0.  The  eigenstructure  of  the  reference  system  is  FX  =  A' A  where  A 
is  a  diagonal  matrix  of  the  eigenvalues  of  F,  A (F)  and  X  is  a  matrix  of  the  associated 
eigenvectors. 
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One  common  design  approach  is  LQ  optimization.  It  has  the  desirable  properties  of 
producing  controllers  which  axe  guaranteed  to  be  stable  through  the  solution  of  a  matrix 
equation.  In  particular,  the  stabilizing  controller  which  minimizes  ||G(s)||2  is  given  by 

tt  =  K2x,  K2  =  -BTM2  (6.6) 

where  M2  >  0  is  the  solution  of  the  algebraic  Riccati  equation 

ATM2  +  M2 A  -  M2BBTM2  +  HtH  =  0.  (6.7) 

For  details,  see  [12]. 

A  stabilizing  controller  which  guarantees  (|G(s)||oo  <  7  is  given  by 

u  —  KooX,  Koo  =  -BT  Mao  (6.8) 

providing  there  exists  Mao  >  0  which  satisfies  the  algebraic  Riccati  equation 

ATMao  +  Mao  A  -  MooBBTMao  +  \MooGGtMoo  +  HtH  =  0.  (6.9) 

72 

For  details,  see  [5]. 

6.2  Static  Controllers 

In  this  section,  the  static  controller  will  be  determined  which  retains  the  r  reference 
eigenvalues  Ar  and  associated  eigenvectors  Xr  where  r  is  the  number  of  measured  outputs. 

Theorem  13  If  Ar  is  observable  from  C,  then  the  unique  static  output  feedback  controller 
which  retains  the  r  reference  eigenvalues  is  given  by 

C(s)  =  Dc  (6.10) 

where 

Dc  =  K0N0  (6.11) 

N0  =  XT{CXTyl  (6.12) 

and  Xr  spans  the  eigenspace  associated  with  the  r  reference  eigenvalues  Ar. 
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Proof:  Since  the  feedback  gain  Dc  retains  [Ar,Xr] 

AcXr  =  (A  +  BDcC)Xr  =  *rAr.  (6.13) 

Also,  note  that 

FXr  =  (A  +  BK0)Xr  =  XrAr.  (6.14) 

Subtracting  the  two  equations  yields 

BDcCXt  =  BK0Xr.  (6.15) 

Since  Ar  is  observable  from  C,  CXr  is  invertible.  Thus, 

Dc  =  K0Xr(CXr)-\  (6.16) 

■ 

The  eigenvalues  of  the  closed  loop  system  consist  of  the  retained  eigenvalues  and  the 
residual  eigenvalues.  Note  that  the  eigenvectors  are  also  retained. 

Theorem  14  The  eigenvalues  of  the  system  are 


Ac  =  Ar  U  A(Ar) 

where 

Ar  =  YT(In  -  N0C)AY 
and  Y  €  7£n*("-r>  is  such  that  CY  =  0  and  YTY  =  7n_r. 

Proof:  Let  T  be  given  by 


T=  Xr  Y 


r 


V 


(CXT)~lC 
YT(In  -  N0C ) 


(6.17) 


(6.18) 


(6.19) 


Note  that  U  and  V  exist  provided  that  CXr  is  invertible  which  is  guaranteed  by  the 
observability  of  Ar.  Thus  T  is  invertible  since  U  and  V  exist. 


T~XACT  = 


A  r  * 

0  Ar 


Ar  =  YT(In  -  N0C)AY. 


(6.20) 

(6.21) 
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Gr(s) 


GP(s) 


Figure  6.1:  Error  system. 


The  following  theorem  forms  the  transfer  function  of  the  system  represented  by  the 
difference  between  the  reference  and  the  projective  systems.  This  system  is  illustrated  in 
Figure  6.1. 

Theorem  15  The  transfer  function  of  the  error  system  E(3 )  =  Gr(s)  —  Gp(s )  can  be 
expressed  as 


E(s)  = 


F 

B 

H  +  EK0 

E 

at 

Gr 

K0Y 

0 

where  Gr  =  YT(In  -  N0C)G. 

Proof:  By  definition,  E(s)  is  given  by 


F 

0 

G 

E(s)  = 

0 

Ac 

-G 

H  +  EK0 

H  +  EDcC 

0 

J 

Applying  the  state  space  transformation, 


T0  = 


E(s )  is  given  by 


E(s)  = 


,  -Xr 

XT 

r  ■  > 

T-\  _ 

0 

In  In 

0  u 

0  V 

F 

0 

-BK0Y 

0 

0 

Ar 

UAY 

-UG 

0 

0 

Ar 

-VG 

(6  22) 


(6.23) 


(6.24) 


(6.25) 


H  +  EK0  0  -EK0Y 


0 
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Removing  the  unobservable  states  yields 


F 

-BK0Y 

0 

E(s)  = 

0 

AT 

-VG 

H  +  EK0 

-EK0Y 

0 

F 

B 

Ar 

VG 

H  +  EK0 

E 

K0Y 

0 

(6.26) 


6.3  Strictly  Proper  Controllers 


In  this  section,  the  use  of  strictly  proper  controllers  for  projective  controls  is  investi¬ 
gated.  It  is  shown  that  the  use  of  strictly  proper  controllers  introduces  free  parameters 
which  can  be  used  to  shape  the  residual  dynamics. 

The  following  theorem  constructs  the  family  of  p-th  order  strictly  proper  controllers 
which  retain  p  reference  eigenvalues  and  eigenvectors. 


Theorem  16  The  set  of  p-th  order  strictly  proper  controllers  which  retain  p  reference 
eigenvalues  and  eigenvectors  is  given  by 


C(s)  = 


Ac 

Sc' 

Cc 

0 

(6.27) 


where 


Ac  —  Ap  —  P0CXp 
Bc  =  P0 
Cc  —  XaXp, 


(6.28) 

(6.29) 

(6.30) 


P0  G  TV>XT  is  a  free  parameter  matrix  and  Xp  spans  the  eigenspace  associated  with  the  p 
reference  eigenvalues,  Ap. 

Proof: 

A  BCe 
BeC  Ac 

ACWP  +  BeCXp  =  WpAp 


_ 1 

* 

Wp 

WP\ 

(6.31) 


Ac  =  WpApW~x  -  BcCXpW -1 


(6.32) 

(6.33) 
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Ac  =  Wp{ AP  -  W;lBcCXp)Wp-\  (6.34) 

Define  PQ  =  W~lBc 

Ac  =  WP(AP  -  P0CXp)W~l  (6.35) 

Bc  =  WpP0.  (6.36) 

BCcWp  +  AXp  =  XpAp  (6.37) 

BCcWp  +  AXp  =  (A  +  BK0)XP  (6.38) 

CcWp  =  KoXp  (6.39) 

Cc  =  K0XpW;\  (6.40) 


Note  that  Wp  represents  a  state  space  transformation  of  the  controller  and  thus  Wp  is 
arbitrary.  Choose  Wv  =  /p,  then 


Ac  =  Ap  -  P0CXp  (6.41) 

BC=P0  (6.42) 

Cc  =  tf0Xp.  (6.43) 


Theorem  17  The  eigenvalues  of  the  closed  loop  system  are 

Ac  =  APU  A(Ar) 


where 

Ar  =  A-  XpPoC. 


Proof: 


Let 


A  BCc 
BcC  Ac 


T  = 


’  X, 

In 

,  f'1  = 

o 

h 

. 

0 

« 

[/» 

~XP 

I 

(6.44) 


(6.45) 


(6.46) 


(6.47) 
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f~lAct 


Ap  * 

0  Ar 


Ar  =  A-  XpP0C 


(6.48) 

(6.49) 


Theorem  18 


The  transfer  function  of  the  error  system  E(s )  =  Gr(s)  —  Gp(s)  is  given  by 


E(s)  = 


F 

B 

H  +  EK0 

E 

Ar 

G 

K0 

0 

Proof: 


E(s)  = 


F 

0 

0 

-G  " 

0 

A 

BCc 

G 

0 

BcC 

Ac 

0 

H  +  EK0 

H 

ECc 

0 

Applying  the  state  space  transformation 


T  = 

O  - 

- 1 

O 

-*p 

Xp 

-In 

In 

II 

In 

0 

In 

0 

t 

o 

0 

h 

0 

0 

In 

v  1 
Xp  j 

yields 


E(s)  = 


F 

0 

-BK0 

0 

0 

Ap 

PoC 

0 

0 

0 

Ar 

G 

H  +  EK0 

0 

-EK0 

0 

Removing  the  unobservable  states  yields 


E(s)  = 


F 

0 

-BK0 

A 

0 

G 

F 

B 

• 

At 

H  +  EK0 

r7» 

H  +  EK0 

-EK0 

0 

Cj 

• 

Ar 


(6.50) 


(6.51) 


(6.52) 


(6.53) 


(6.54) 
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6.4  Proper  Controllers 

In  this  section,  the  use  of  proper  controllers  (that  is,  dynamic  cr-;itrollers)  for  projective 
controls  is  investigated.  It  is  shown  that  the  use  of  proper  controllers  allows  more  of  the 
reference  dynamics  to  be  retained  and  introduces  free  parameters  which  can  be  used  to 
shape  the  residual  dynamics. 

The  following  theorem  constructs  the  family  of  p-th  order  proper  controllers  which 
retain  r  +  p  reference  eigenvalues  and  eigenvectors. 

Theorem  19  The  set  of  all  p-th  order  proper  controllers  which  retain  r  +  p  reference 
eigenvalues  and  eigenvectors  is  given  by 


where 


(6.55) 


C(s)  = 


'  Ac 

Bc  ' 

Cc 

Dc 

Ac  =  Ap  +  P0CFB0  (6.56) 

Bc  =  P0CF(N0  -  B0P0 )  -  A pP0  (6.57) 

Cc  =  K0B0  (6.58) 

Dc  =  K0(N0  -  B0P0 ),  (6.59) 


B0  =  (I  —  N0C)XP,  P0  €  RpXT  is  a  free  parameter  matrix  and  XT  a.nd  Xp  span  the  reference 
eigenspaces  associated  with  the  reference  eigenvalues,  Ar  and  Ap. 

Proof: 


AcXr  =  XrAr 


A  4-  BDcC 

BCc 

’  xP 

Xp  Xr 

Ap  0 

BcC 

Ac 

Wp 

wr 

£ 

1 _ 

0  Ar 

(6.60) 

(6.61) 


Ae  Wp  +  BcCXp  =  WP  Ap  ( 6. 62 ) 

AcWr  +  BcCXr  =  WrAr  (6.63) 

BCcWp  +  (.4  +  BDcC)Xp  =  .YpAp  =  ( A  +  BI<0  )XP  (6.64) 

BCcWr  +  {A  +  BDcC)Xr  =  .YrAr  =  (.4  +  BK0)Xr  (6.65) 
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Ac  Bc  WB  Wr 


Cc  Dc  CXV  CXr 


WpA.p  Wr  Ar 
KoXp  KoXr 


Define  L  =  W~xWr 


Ac  Bc 
Cc  Dc 


Ac  Bc 
Cc  Dc 


Ac  Be 

Cc  Dc 


Wp  0 


WpAp  Wr  Ar  Wp  Wr 


K0Xp  KaXr  I  CXp  CXr 


(6.66) 


(6.67) 


WpAp  WpLAr  Wp  WpL 


K0XP  K0Xr  CXp  CXr 


Av  LAr 


0  Im  I  KoXp  K0Xr  I  CXp  CXr 


Ip  L 


V  0 

0  Ir 


(6.68) 


(6.69) 


Note  that  Wp  represents  a  state  space  transformation  of  the  controller  and  thus,  Wp  is 
arbitrary. 


h  L 

CXp  CXr 


Ip  +  LA-'CXp  -LA-1 
-A-'CX,  A-1 


A  =  CXr  -  CXpL. 


Define  PQ  =  L{CXr  -  CXPL)~\  then 


(6.70) 


(6.71) 


h  L 

CXp  CXr 


Let  Wp  =  (Ip  +  P0CXp) 


A-1  =  (CXr)~1(IT  +  CXpPo). 


Ip  +  P0CXp 


-(CXr)~lCXp(Ip  +  PoCXp)  (CXr)~\Ir  +  CXpP0) 


Ac  Bc 
Cc  Dc 


Ac  Bc 
Cc  Dc 


(Ip  +  P0CXp)Ap  PoCXr  Ar 


K0Xr 


(6.72) 


(6.73) 


KoXp  KoXr  I  I  -(CXr)~l  CXp  (CXr)-l(Ir  +  CXpPo) 

(6.74) 

Ap  +  P0CF(XV  -  N0CXp)  P0CF(N0CXp  -  Xp)P0  +  P0CFN0  -  A pP0 
K0(XP  -  N0CXp)  K0(N0CXp  -  Xp)P0  +  K0N0 

(6.75) 


Ac  Be 

Cc  Dc 


Ap  +  PpCFBo  PoCF(N0  -  B0P0)  -  A VP0 
KoB0  K0(N0  —  B0P0) 


(6.76) 
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Theorem  20  The  eigenvalues  of  the  closed  loop  system  are 

Ac  =  ApUArUA(ir) 

where 

Ar  =  Ar  +  B0P0CAY 

and  Y  6  TZnx^n~T^  is  such  that  CY  =  0  and  YTY  =  J„_r. 
Proof: 


Ac  = 


A  +  BDcC  BCc 


(6.77) 


(6.78) 


(6.79) 


-  =  Xp  Xr  Y 

[lp  +  P0CXV  P0CXT  0 

.  n  -PoC  iP 

u 

r_1  =  .  =  ( cxr)-lc(in  +  xpp0c )  -(cxr)-lcxp 

L  J  Yt  +  Yt(B0P0-N0)C  - YTBo 

Ap  0  * 

f-lAcf  =  0  Ar  * 

0  0  Ar 

Ar  =  YT(In  +  ( B0P0  -  N0)C)AY  =  Ar  +  YtB0P0CAY. 


(6.80) 


(6.81) 


(6.82) 

(6.83) 


Theorem  21  The  transfer  function  of  the  error  system  E(s)  =  Gr{s)  —  Gp(s)  is  given  by 

[  F  I  B  1  F  Ar  Gr  +  YT B0P0CG  1 


E(s)  = 


H  +  EK0  E  K0Y 


(6.84) 


Proof: 


E(s)  = 


F 

0 

0 

-G  ' 

0 

A  +  BDcC 

BCc 

G 

0 

BcC 

Ac 

0 

H  +  EK0 

H  +  EDcC 

ECc 

0 

(6.85) 
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Applying  the  state  space  transformation 

In  ~Xp 


t  = 


j-1 


-XT  -Y 
0  Xp  Xr  Y 

0  Ip  +  P0CXp  PoCXr  0 

In  In  0 

0  ~PoC  /p 

0  (CXr)~lC(In  +  XpP0C)  ~{CXr)~xCXp 
0  YT  +  YT(B0P0  -  No)C  -YtBo 


yields 


E(s)  = 


F 

0 

0 

-BK0Y 

0 

0 

Ap  0 

* 

* 

0 

0  Ar 

* 

* 

0 

0 

0 

Ar 

Gr  +  YTB0P0CG 

H  +  EK0 

0 

0 

-EK0Y 

0 

Removing  the  unobservable  states  yields 

E(s)  = 

which  is  equivalent  to 

E(s)  = 


F 

-BK0Y 

0 

0 

Ar 

Gr  +  YrB0P0CG 

H  +  EK0 

-EK0Y 

0 

F 

B 

H  +  EK0 

E 

Ar 

Gr  +  YTB0P0CG 

K0Y 

0 

(6.86) 


(6.87) 


(6.88) 


(6.89) 


(6.90) 


6.5  The  Auxiliary  Minimization  Problem 

The  previous  sections  established  projective  controllers  in  terms  of  free  parameters. 
These  free  parameters  affected  the  residual  dynamics  of  the  system  and  thus  are  an  impor¬ 
tant  part  of  the  design  problem.  However,  no  straightforward  method  exists  for  choosing 
these  parameters.  In  this  section,  an  auxiliary  minimization  problem  is  formulated  in  order 
to  determine  the  free  parameters  of  the  dynamic  projective  controllers. 
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The  auxiliary  minimization  problem  is  as  follows:  Determine  P0  such  that  ||P(s)||f  is 
minimized. 

Thus,  it  is  natural  to  apply  the  results  of  Chapter  4  to  determine  the  optimal  parameters 
P0.  In  the  dynamic  controller  cases,  E(s)  is  given  by 

E(s)  =  El(s)E2(s]P0).  (6.91) 


Note  that  E(s)  is  a  system  of  order  2 n.  Thus,  the  auxiliary  minimization  problem  is  twice 
the  order  of  the  original  problem.  Also,  the  structure  of  E(s)  as  shown  in  (6.91)  cannot 
be  exploited  to  simplify  the  problem.  These  problems  do  not  prevent  the  FH  optimization 
method  from  being  solved  and  necessary  conditions  for  an  optimal  solution  can  readily  be 
found. 

However,  since  Ei(s )  is  independent  of  Pa,  it  acts  only  as  a  weighting  factor  and  does  not 
influence  the  stability  or  the  robustness  properties  of  the  system.  Thus,  as  a  computational 
simplification,  the  auxiliary  minimization  problem  is  reduced  to 


min||£2(s)||F. 

*o 


(6.92) 


Note  that  this  is  not  the  same  as  minimizing  E(s )  and  in  general  will  give  a  different 
solution.  However,  it  still  represents  minimizing  a  weighted  sum  of  the  error  system  E(s). 
In  the  case  of  strictly  proper  projective  controllers,  P2(s)  is  given  by 


A  -  XVP0C 

G  " 

Ar 

Br 

K° 

0 

Gr 

0 

(6.93) 


Thus,  the  auxiliary  minimization  problem  is  to  minimize  over  P0, 


J(Po)  —  Tr  PQ, 


(6.94) 


subject  to  the  two  Lyapunov  equations  which  define  the  controllability  grammian  P  and 
the  observability  grammian  Q: 


ArP  +  PAtt  +  BrBj  =  0 

(6.95) 

AjQ  +  QAr  +  C?Cr  =  0. 

(6.96) 
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The  necessary  conditions  for  an  optimal  solution  from  Chapter  4  are  (6.95)-(6.96)  and 

AjL  +  LAr  +  Q  =  0  (6.97) 

ArM  +  MAj  +  P  =  0  (6.98) 

Xp(LP  +  QM)Ct  =  0  (6.99) 

where  M  and  L  are  Lagrange  multipliers  for  the  constraints  (6.95)-(6.96). 

The  gradient  of  J  with  respect  to  P0  can  be  computed  for  arbitrary  P0  as 

=  -2Xj(QL  +  MP)CT  (6.100) 

where  P,  Q,  L  and  M  solve  the  Lyapunov  equations  (6.95)-(6.98). 

Thus,  a  steepest  descent  algorithm  can  be  implemented  by  iteratively  solving 

*"“*-‘(35)  (6101) 

for  the  optimal  Pa. 

This  method  of  determining  free  parameters  is  employed  in  the  design  example  of 
Chapter  8. 


CHAPTER  7 


DECENTRALIZED  PROJECTIVE  CONTROLS 

In  this  chapter,  projective  controls  method  for  the  case  of  decentralized  controllers  is 
examined.  Much  of  the  work  presented  in  this  chapter  is  developed  in  [32,19,18].  As  in 
the  centralized  case,  the  presentation  here  focuses  on  a  transfer  function  point  of  view. 
In  addition,  a  new  formulation  of  the  proper  controller  case  is  developed,  which  avoids 
the  requirement  of  transforming  the  plant  into  a  specific  form,  and  results  for  the  case  of 
strictly  proper  controllers  are  introduced. 

The  main  focus  of  this  chapter  is  to  show  that  the  free  parameters  of  the  decentralized 
projective  controllers  can  be  treated  in  a  manner  similar  to  the  centralized  case.  Thus,  FH 
optimization  is  suitable  for  determining  the  free  parameters  of  the  controller.  However,  in 
order  to  efficiently  apply  the  FH  norm  minimization  approach  to  the  disturbance  rejection 
problem,  a  transformation  developed  in  [32]  will  be  applied  to  the  system  that  reduces  the 
closed  loop  system  to  a  form  which  is  linear  in  the  free  parameters. 

Consider  the  system  with  l  feedback  channels 


The  goal  is  to  determine  the  /  decentralized  controllers  Ci(s)  to  achieve  certain  performance 
and  robustness  goals. 

7.1  The  Reference  System 

Given  the  l  state  feedback  controllers  u,  =  KaiX,  the  reference  system  is  given  by 
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where 

/ 

F  =  A  +  Y,B'Ko'-  (7-3) 

1=1 

7.2  Static  Controllers 


For  completeness,  consider  first  the  use  of  static  output  feedback. 

Theorem  22  The  l  decentralized  static  controller  which  retains  r  reference  eigenvalues  is 
given  by 

Ci(s)  =  Ki  (7.4) 

where 

K,  =  KoiNoi  (7.5) 

Not  =  XriCiX  r)-1  (7.6) 

and  XT  spans  the  eigenspace  associated  with  the  r  reference  eigenvalues  Ar. 


The  transfer  function  of  the  closed  loop  system  is 


where 

( 

Ac  =  A  +  Y,BxKtC,.  (7.10) 

i=i 

Hence,  with  static  controls,  there  is  no  freedom  available  for  the  next  design  phase  within 
the  class  of  projective  controls.  This  solution,  however,  serves  as  a  foundation  for  expanding 
the  admissible  controllers  by  allowing  dynamic  controllers  of  given  order.  The  advantages 
are  twofold:  more  eigenvalues  can  be  retained  and  free  parameters  are  introduced  to  shape 
the  residual  dynamics  to  achieve  disturbance  rejection.  The  proper  controller  and  strictly 
proper  controller  cases  Eire  treated  seperately. 
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7.3  Strictly  Proper  Dynamic  Controllers 

Theorem  23  The  set  of  all  l-channel,  p-th  order  decentralized  strictly  proper  controllers 
which  retain  the  p  reference  eigenvalues  is  given  by 

Ad  Bd 

Ci(s)  =  — —  (7.11) 

ca  0 

where 


(7.11) 


Ad  —  Ap  —  P,CiXp 

(7.12) 

fco 

2. 

II 

(7.13) 

Cd  —  XoiXp , 

(7.14) 

Pi  €  RpXr  are  free  parameter  matrices  and  Xp  spans  the  reference  eigenspaces  associated 
with  the  p  reference  eigenvalues,  Ap. 

Proof:  ■ 

The  transfer  function  of  the  closed  loop  system  is  then 


G(s)  = 


Ac  G 
H  0 


” 

A 

BiK0iXp 

BiK^Xp 

G 

PxC i 

Ap  -  PlC,Xp 

0 

0 

p2c2 

0 

Ap  -  P2C2Xp 

0 

11 

0 

0 

0 

(7.15) 


Note  that  this  system  is  linear  in  the  free  parameters.  The  freedom  in  the  free  pa¬ 
rameters  P\  and  P2  can  be  used  to  shape  the  residual  dynamics  and,  in  the  disturbance 
attenuation  problem,  to  shape  the  system’s  transfer  function.  This  linear  representation 
(7.15)  enables  the  efficient  application  of  the  FH  norm  approach. 

7.4  Proper  Controllers 

Theorem  24  The  set  of  all  l-channel ,  p-th  order,  decentralized  proper  controllers  which 
retain  r  +  p  reference  eigenvalues  is  given  by 

Ad  Bd 

Ct(s)  =  — -  (7.16) 

Cd  o 


(7.16) 


Ad  =  A  p  +  PiCiF  B0i 

(7.17) 

Bd  =  P,CtF(Not  -  Boi Pi)  ~  ApP, 

(7.18) 

Cd  =  KoiBoi 

(7.19) 

Dd  =  KoiiNoi  -  BoiPi), 

(7.20) 

Boi  =  (I  —  NoiCi)Xp,  P{  €  Rpxr  are  free  parameter  matrices  and  XT  and  Xp  span  the 
reference  eigenspaces  associated  with  the  reference  eigenvalues,  Ar  and  Ap. 

The  transfer  function  of  the  closed  loop  system  is 


<?(*)  = 


Ac  G 
H  0 


Ad 

B\Cc\ 

B2Cc2 

BdCi 

Ad 

0 

Bc2C2 

0 

AC2 

H 

0 

0 

where 


Ad  —  Ac  4-  B\Dc\C\  +  B2Dc2C2. 

Transforming  the  closed  loop  system  by  the  similarity  transformation 


/„  0  0 

T=  P,C,  Ip  0 

P2C2  0  Ip 


In  0  0 
T-1  =  —P\C\  Ip  0 

-P2C2  0  Ip 


yields 


Ac  G 


Ac  B\Cci  B2Cc2  G 

P\E\  Ap  +  PtGu  -PG12  -P\C\G 


(7.21) 


(7.22) 


(7.23) 


Ci{FNoiCx 
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Thus,  the  expressions  derived  for  Ac,  G  and  H  ail  exhibit  a  linear  dependence  on  the  free 
parameter  matrices  P\  and  P2.  (H  is  in  fact  independent  of  the  free  parameter  matrices.) 
This  linear  dependence  can  now  be  utilized  to  determine  suitable  P\  and  P2  (and  thus  the 
dynamic  controllers)  to  achieve  disturbance  attenuation  by  minimizing  the  FH  norm. 
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CHAPTER  8 

A  CONTROLLER  DESIGN  EXAMPLE 

The  methods  presented  in  the  previous  chapters  are  now  applied  to  the  problem  of 
controlling  a  flexible  structure.  The  problem  explored  in  this  design  example  is  the  im¬ 
provement  of  the  disturbance  attenuation  properties  of  a  flexible  structure  using  low  order, 
robust  control. 

The  structure  considered  here  is  the  cruciform-shaped  structure  shown  in  Figure  8.1. 
Note  that  this  is  a  realistic  design  problem  since  the  plant  model  used  in  this  design  problem 
is  taken  from  an  actual  structure  and  the  model  has  been  verified  experimentally  [10].  In 
addition,  the  structure  exhibits  many  properties  typical  of  flexible  structures  including 
lightly  damped  modes,  decentralized  structure  and  multiple  controls  and  measurements. 
This  structure  has  also  been  studied  in  [28]  and  [26]. 

The  procedure  used  to  design  the  controller  is  as  follows.  First,  after  determining  that 
the  system  exhibits  weak  coupling  between  its  x-axis  and  y- axis  dynamics,  the  problem  is 
separated  into  two  parts:  the  design  of  the  x-axis  controller  and  the  design  of  the  y-axis 
controller.  Second,  a  reference  state  feedback  controller  is  determined  which  achieves  the 
desired  disturbance  attenuation.  Third,  by  selecting  the  proper  modes  to  retain,  the  strictly 
proper  projective  controller  is  determined.  Fourth,  the  free  parameters  in  the  projective 
controller  are  determined  by  applying  FH  optimization  to  solve  the  auxiliary  minimization 
problem.  Finally,  the  properties  of  the  resulting  closed  loop  system  are  shown. 

8.1  Description  of  the  Structure 

The  structure  considered  is  a  45-foot  lattice-type,  lightweight  (5  lb),  flexible  beam  with 
fixed  base  and  free  tip  shown  in  Figure  8.1  [10]. 

The  control  u  (E  TV  consists  of  torques  applied  at  the  base  of  the  structure  about  the 
x  and  y  axes.  The  disturbance  w  €  TV  is  generated  by  an  x-y  translation  applied  to  the 
base  where  the  z-axis  is  taken  to  be  the  axis  of  the  cruciform.  Measurements  of  the  system 
y  €  are  obtained  from  an  x-y  axis  gyro  and  accelerometer  sensors  located  at  the  tip 
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Figure  8.1:  The  cruciform  structure. 

and  base  of  the  structure.  The  controlled  output  z  6  TZ4  is  the  position  measurement  at 
the  tip  and  base  of  the  structure. 

The  system  model  was  obtained  by  truncating  the  lowest  20  modes  of  the  finite  element 
model.  This  results  in  a  fortieth  order  state  space  model.  The  model  parameters  can  be 
found  in  Appendix  D  [10]. 

The  cruciform  structure  was  observed  to  exhibit  significant  decoupling  between  the  x 
and  y  axis  dynamics.  Thus,  the  design  problem  has  been  split  into  two  parts:  The  x-axis 
dynamics  and  the  y-axis  dynamics.  The  inputs  and  outputs  of  the  structure  are  divided 
as  shown  in  Table  D.l  in  Appendix  D.  For  more  details  see  [26]. 

The  x-axis  plant  model  is  obtained  by  applying  a  balanced  model  reduction  on  the 
original  system  using  only  the  x-axis  inputs  and  outputs.  The  eigenvalues  of  the  resulting 
twelth  order  model  are  shown  in  Table  8.1  and  show  typical  flexible  structure  properties, 
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Table  8.1:  Modes  of  the  Cruciform  Model. 


Real 

Imaginary 

Frequency 

Damping 

-4.2471e-03 

±8.4941e-01 

8.  '«542e-01 

5.0000e-03 

-3.4723e-02 

±7.0339e+00 

7.0340e+00 

4.9365e-03 

-3.0800e-02 

±7.7296e+00 

7.7297e+00 

3.9846e-03 

-1.9863e-01 

±1.0458e+01 

]  .0460e+01 

1.8989e-02 

-5.1863e-01 

±2.5925e+01 

2.5930e+01 

2.0001e-02 

-8.2976e-01 

±4.6092e+01 

4.6100e+01 

1.7999e-02 

i.e.,  lightly  damped  and  closely  packed  low  frequency  modes.  The  y-axis  dynamics  were 
treated  similarly  but  are  not  shown  here. 

Note  that  model  reduction  was  used  only  to  remove  very  weakly  controllable  and  ob¬ 
servable  modes  in  order  to  avoid  neglecting  modes  which  may  be  important  in  the  design 
of  the  controller.  This  is  possible  since  the  projective  controls  method  allows  the  use  of  a 
high  order  model  without  requiring  a  high  order  controller. 

o.2  Design  of  the  Controller 

The  goal  of  this  section  is  to  design  a  controller  which  improves  the  disturbance  atten¬ 
uation  of  the  structure.  The  measure  of  disturbance  attenuation  used  here  will  be  the  Hoo 
norm  since  it  represents  the  maximum  gain  of  the  system  from  disturbance  to  controlled 
output.  Thus,  a  system  such  that  'Hoc  <  7  is  said  to  have  disturbance  attenuation  7. 

From  the  frequency  response  shown  in  Figure  8.2,  it  can  be  seen  that  the  open  loop 
system  has  a  disturbance  attenuation  of  7  =  -20  dB.  In  this  example,  the  design  goal  will 
be  to  improve  the  disturbance  attenuation  of  the  structure  to  7=-40  dB  using  a  low  order, 
robust  controller.  This  will  be  done  using  the  projective  controls  method. 

The  first  step  in  designing  a  projective  controller  is  to  find  a  state  feedback  which  solves 
the  disturbance  attenu  ‘-on  problem.  This  forms  the  reference  system  for  the  problem. 
The  following  theorem  1  [26]  yields  the  desired  controller. 


Theorem  25  Given  a  plant 


A 

G 

B 

P(s)  = 

H 

0 

0 

In 

0 

0 

(8.1) 


where  there  exists  an  a  such  that  G  =  Ba  and  a  state  feedback  controller  C(s )  =  K  where 


K  =  -R~lBTP 


(8.2) 


and  P  satisfies  the  Riccati  equation 

AtP  +  PA-  PBtR~xBP  +  HtH  =  0 


(8.3) 


where  R  satisfies 

aTRa  <  72,  (8.4) 

then  the  closed  loop  system  satisfies  the  bound  ||<!t(.s)||oo  5:  7- 

Making  the  approximation  D  «  Ba  where  a=0.56,  Theorem  25  can  now  be  applied 
to  yield  the  state  feedback  controller.  The  eigenvalues  of  the  resulting  closed  loop  system, 
A (F),  are  given  in  Table  8.2.  The  frequency  response  of  the  reference  system  is  shown  in 
Figure  8.3  and  confirms  that  the  disturbance  attenuation  of  the  reference  system  is  7  = 
-40  dB.  The  robustness  margins  are  [0,  00)  for  gain  and  ±90°  for  phase.  Thus,  this  state 
feedback  forms  an  acceptable  reference  system  for  the  projective  controls  method. 

First,  a  static  projective  controller  is  considered.  Since  there  are  r  =  4  outputs,  two 
complex  pair  modes  can  be  retained.  However,  experimentation  with  possible  retained 
modes  does  not  yield  any  closed  loop  systems  with  the  required  disturbance  attenuation. 
Since  no  design  freedom  exists  for  the  static  projective  case,  nothing  further  can  be  done 
in  this  case. 

Since  the  static  projective  controls  approach  was  inadequate,  the  strictly  proper  projec¬ 
tive  controls  approach  is  now  considered.  Using  a  second  order  ( p  =  2)  controller  (higher 
order  controllers  can  be  considered  later  if  this  is  inadequate),  one  complex  reference  mode 
pair  can  be  retained.  Referring  to  Figure  8.3,  it  is  apparent  that  the  lowest  frequency 
mode  of  the  reference  system  is  the  most  important  in  retaining  disturbance  attenuation. 
This  is  the  mode  labeled  “a’’  in  Table  8.2.  Thus,  choose  the  retained  mode  to  be  AP  =  {a}. 


frequency  (radians/second) 


Figure  8.2:  Frequency  response  of  the  open  loop  system. 


Table  8.2:  Modes  of  the  Reference  System 


Real 

Imaginary 

Frequency 

Damping 

a 

-1.0856e-01 

±8.5272e-01 

8.5960e-01 

1.2629e-01 

b 

-3.4723e-02 

±7.0339e+00 

7.0340e+00 

4.9365e-03 

c 

-3.0801e-02 

±7.7296e+00 

7.7297e-f00 

3. 984 8e- 03 

d 

-1.9902e-01 

±1.0458e+0I 

1.0460e+01 

1.9026e-02 

e 

-5.1870e-01 

±2.5925e+01 

2.5930e+01 

2.0004e-02 

f 

-8.2977e-01 

±4.6092e+01 

4.6100e+01 

1.7999e-02 
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The  choice  of  Ap  determines  the  family  of  strictly  proper  projective  controllers  which 
retain  Ap  parameterized  by  PQ.  To  select  Pa,  the  approach  of  Section  6.5  is  applied.  Using 
a  gradient  method,  the  Frobenius-Hankel  norm  of  £2(5)  was  minimized.  The  resulting 
controller  is  given  in  Appendix  E.  The  frequency  response  of  the  dynamic  controller 
associated  with  this  choice  of  Pa  is  given  in  Figure  8.4. 

8.3  Evaluation  of  Design 

The  final  step  is  to  evaluate  the  design  by  examining  the  properties  of  the  closed  loop 
system.  The  eigenvalues  of  the  closed  loop  system  are  shown  in  Table  8.3.  Note  that 
Ap  =  {a}  has  been  retained.  The  disturbance  attenuation  of  the  full  system  is  7  =  -40  dB 
as  seen  in  Figure  8.5.  Thus,  the  disturbance  attenuation  goal  of  the  design  has  been  met 
using  a  second  order  controller.  The  stability  margins  of  the  closed  loop  system  are  [0,  40 
dB]  in  gain  and  ±70°  in  phase.  While  these  stability  margins  are  smaller  than  those  of  the 
reference  system,  they  are  still  quite  large.  However,  if  these  margins  are  not  satisfactory, 
a  higher  order  controller  may  be  considered. 

To  demonstrate  the  disturbance  attenuation  achieved  by  this  design  in  the  time  domain, 
the  time  response  of  the  system  to  a  disturbance  impulse  is  computed.  The  open  loop 
response  is  given  in  Figure  8.6.  Note  the  low  damping  of  the  low  frequency  mode.  For  the 
system  controlled  by  the  design  given  above,  the  response  is  given  in  Figure  8.7.  In  this 
case,  the  damping  on  the  low  frequency  mode  has  increased  dramatically. 

A  second  order  controller  was  also  designed  for  the  y-axis  dynamics  of  the  system  in  a 
similar  manner.  The  resulting  closed  loop  system  with  decentralized  controls  was  seen  to 
be  stable  and  retain  the  desired  disturbance  attenuation  properties. 

Thus,  the  disturbance  attenuation  of  a  flexible  system  modeled  with  a  fortieth-order 
model  has  been  improved  significantly  using  two  second  order,  decentralized  controllers. 
The  use  of  two  second  order  controllers  is  quite  an  improvement  over  the  full  fortieth  order 
controller.  Thus,  the  projective  controls  method  with  FH  optimization  has  proven  to  be 
an  effective  design  method  in  a  realistic  problem  setting. 
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Table  8.3:  Modes  of  the  Closed  Loop  System 


Real 

Imaginary 

Frequency 

Damping 

-2.5695e-01 

±7.7933e-01 

8.2059e-01 

3.1313e-01 

-1.0856e-01 

±8.5272e-01 

8.5960e-01 

1.2629e-01 

-3.4630e-02 

±7.0333e+00 

7.0334e+00 

4.9237e-03 

-3.0564e-02 

±7.7291e+00 

7.7292e+00 

3.9543e-03 

-1.8496e-01 

±1.0423e+01 

1.0425e+01 

1.7742e-02 

-5.3975e-01 

±2.5960e+01 

2.5965e+01 

2.0787e-02 

-8.1569e-01 

±4.6082e+01 

4.6089e+01 

1.7698e-02 

Figure  8.5:  Frequency  response  of  the  closed  loop  system. 
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CHAPTER  9 
CONCLUSIONS 

This  work  has  developed  methods  of  designing  low  order,  decentralized  controllers 
which  give  robust  performance  in  the  closed  loop  system. 

The  robustness  properties  of  the  /H00  norm  make  it  an  attractive  choice  for  computing 
low  order  optimal  controllers.  However,  no  methods  previously  existed  to  solve  such  a 
problem.  Necessary  conditions  for  the  solution  of  this  problem  shown  in  this  work,  however, 
still  do  not  lead  to  effective  methods  for  the  computation  of  such  controllers. 

The  FH  norm,  however,  was  shown  to  be  a  suitable  norm  for  use  in  the  design  of  low 
order  controllers  based  on  its  ability  to  quantify  both  robustness  and  performance  issues. 
In  addition,  it  poses  a  much  simpler  computational  problem  which  makes  it  suitable  for 
application  to  complex  problems. 

The  FH  optimization  method  was  used  in  the  design  of  low  order  controllers  in  two 
distinct  ways.  First,  necessary  conditions  for  FH  optimal  controllers  were  derived.  This 
led  to  an  optimization  algorithm  which  was  demonstrated  through  an  example.  In  the 
second  approach,  projective  controls  were  used  in  the  first  phase  of  the  design.  However, 
the  projective  controllers  contain  free  parameters  which  influence  the  residual  dynamics  of 
the  closed  loop  system.  FH  optimization  was  then  used  to  choose  the  free  parameters  to 
match  the  closed  loop  system  as  closely  as  possible  with  the  reference  system. 

As  an  example  of  this  approach,  a  design  example  was  shown  which  used  projective 
controls  and  FH  optimization  to  determine  the  control  of  a  flexible  structure. 

With  the  recent  development  of  new  methods  to  solve  the  full  order  Hoo  problem,  the 
problem  of  determining  full  order,  centralized  controllers  can  be  considered  to  be  close  to  a 
solved  problem.  However,  the  situation  is  far  different  for  the  case  of  low  order  controllers. 
The  problem  of  determining  the  best  low  order  controller  is  still  an  open  problem.  A 
variety  of  techniques  are  necessary  to  adequately  design  such  controllers.  This  work  has 
shown  some  new  ways  to  attack  this  problem. 


64 


APPENDIX  A 

PROPERTIES  OF  THE  TRACE  FUNCTION 


Definition  7  Given  A  E  'JZnXn 

Tr  A  =  £>« 

1=1 

In  the  following,  assume  A  E  7Znxn,  B  E  'Rnxm,  C  E  7lmxn. 


Property  1 


Tr  BC  =  Tr  CB. 


(A.l) 


(A.2) 


Property  2 


Tr  A  =  Tr  AT . 


(A. 3) 


Property  3 


Property  4 


TrX  =  X>i(A). 

1  =  1 


Tr  BtC  =  C. 


(A. 4) 


(A. 5) 
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APPENDIX  B 

FH  OPTIMAL  STATE  FEEDBACK  CONTROLLERS 

This  theorem  shows  that  the  unconstrained  FH  optimal  controller  is  not  equivalent  to 
a  static  state  feedback  controller  as  is  the  case  for  an  LQ  controller. 

Theorem  26 

inf  G(s)  <  inf  G(s).  (B.l) 

u=  Kx 

Proof:  The  minimum  principle  is  as  follows:  Let 

x  =  f(x,u,t),  x(t0)  =  x0.  (B.2) 


Find  the  control  u  which  minimizes  the  criterion 

rtf 

J=  !(x,u,t)dt +  m(x(tf)). 
Jto 

Define  the  Hamiltonian  of  the  system  as 

H  =  l  +  pTf. 

Then  the  necessary  conditions  for  a  minimum  are 

x  =  /(x,u,t),  x(t0)  =  x0 
p=-S7*H,  p(tj)  =  S7xjm 
0  =  S7uH. 

Now  consider  the  LTI  system, 

x  —  Ax  -f  B2u,  x(t0)  =  x0 
z  =  Cjx  +  D\2w 
y  =  C2x  +  D2\u. 

For  the  FH  norm,  define  the  following  criterion 

J  =  -  [  tzT zdt. 

2  Jt0 


(B.3) 


( B  .4 ) 


( B  .5 ) 


(B.6) 

(B.7) 

(B.S) 


(B.9) 
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Then, 

H  =  ^t(xTCfCix  +  uTu)  +  pr(Ax  +  B2u). 

This  yields 

x  =  Ax  +  B2u,  x(tQ)  =  x0 

p  =  —tC^C\x  —  ATp,  p(tj )  =  0 
0  =  tu  +  Bjp. 

Thus,  the  optimal  control  is 

it  =  —t~lB2p. 

To  convert  the  control  to  feedback  format,  assume  p(t)  has  the  form, 

p(t)  =  tP(t)x(t). 

Then, 

p  =  Px  +  tP  +  tPx. 

Substitution  yields 

-  P  =  AtP  +  PA  +  CfCj  -  PB2B\ P  +  l-p. 
Thus,  the  optimal  control  is  given  by 

u  =  K(t)  x 


where 

K{t)  =  -BjP(t). 

Note  that  as  t0  —*  — oo,  (B.15)  does  not  have  a  steady-state  solution.  To  see  this, 
that  a  steady-state  solution  P  exists.  Then, 

P  =  0^-P  =  C^P  =  Ct 

t 

which  is  a  contradiction.  Thus,  a  steady-state  solution  does  not  exist. 

This  shows  that  the  unconstrained  FH  optimal  controller  is  in  fact  a  time 
controller  and  not  equivalent  to  a  static  state  feedback  controller. 


(B.10) 

(B.ll) 

(B.12) 

(B.13) 

(B.  14) 

(B.15) 

(B.16) 

(B.17) 

assume 

(B.1S) 

varying 
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APPENDIX  C 

FH  OPTIMIZATION  ALGORITHM 


The  following  is  a  Matlab  macro  which  implements  the  algorithm  for  determining  the 
FH  optimal  controller. 

'/.  FH  minimization 

*/. 

[nn,nz]  =  size(a) ;  [nz,nm] 

[nr,nz]  =  size(c) ;  [nz,nq] 

[ns,nz]  =  size(h) ;  [np,nz] 

'/(Create  extended  system 

A  *  [a  zeros (nn, np) ; zeros (np,nn)  zeros(np)]; 

B  *  [b  zeros (nn,np) ;  zeros(np.nm)  eye(np)] ; 

C  *  [c  zeros(nr,np) ;  zeros(np,nn)  eye(np)] ; 

G  =  [g;  zeros (np, nq) ] ; 

H  *  [h  zeros(ns ,np)] ; 

D  *  [d;  zeros (np, nq) ] ; 

E  =  [e  zeros(ns ,np)] ; 

K  *  [Dc  Cc;Bc  Ac] ; 

tol  *  le-6;  err  =  1;  ep  =  .1; 
yl  *  □  ;  y2  =  □  ; 
while  err>tol. 


Ae  = 

A 

+ 

B*K*C 

Be  = 

G 

+ 

B*K*D 

Ce  = 

H 

+ 

E*K*C 

if  any(real (eig(Ae) )>0) , 

disp(’The  system  is  unstable’), 

end 

'/.Calculate  P  and  Q 
P  =  lyap(Ae , Be*Be ’ ) ; 

Q  =  lyap(Ae’ ,Ce’*Ce) ; 

J  =  trace(P*Q) ,  yl  =  [yl  J] ; 

'/.Calculate  L  and  M 
L  =  lyap(Ae’ ,q) ; 

M  ■  lyap(Ae.P) ; 


=  size(b); 

=  size(g) ; 

=  size(Ac); 


dK  *  B’*(L*P+q*M)*C’  +  B ’ *L* (G+B*K*D) *D ’  +  E ’ * (H+E*K*C) *M*C ’ ; 
z  *  . l*norm(K)/norm(dK) ; 
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if  z  <  1,  ep  =  z; 
else,  ep  =  1; 
end 

K  *  K  -  ep*dK ; 

err  *  norm(dK) ,  y2  —  [y2  err] , 

end 
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APPENDIX  D 
MODEL  PARAMETERS 


The  state  space  model  which  describes  the  cruciform  structure  of  Figure  8.1  is  given  in 


[10]; 


x  =  Ax  -f  Bu  4-  Gw 
y  =  Cx-\-Fu-\-Dw 
z  —  Hx 


where, 


— 2£f2 

-ft2  j 

B  = 

'  bl  ' 

I 

o 

0 

CT  0 

-2 ca&  -c0ft2 


F 


0 

CaBL 


D  = 


0 

CaDL 


H  = 


0  Ca 


ft  =  diag[  .00025 

.00001 

.00083 

.84942 

.87782 

2.2420 

7.0290  7.2100 

7.3300 

7.6228 

7.7327 

7.9790 

8.0682 

8.5011 

10.460 

19.910  20.3S0 

25.930 

42.760 

46.100 

£  =  diag[  .000  .000  .000  .005  .020  .005  .005  .005  .005  .005 

.005  .020  .005  .016  .019  .020  .005  .020  .011  .018  ]. 


7.144e  -  07 

8.829e  -  13 

1.417c -09 

3.503e  -  02 

— 6.267e  -  14 

—  2.170e  -  08 

3.807e  -  02 

— 6.82  le  -  11 

—  1.651c  -  11 

1  380e  -  1 1 

— 3.549e  -  07 

—  1 .698e  -  08 

-1.564c  -  02 

2.570e  -  07 

—  1.522e  -  07 

— 6.110e  -  03 

3.654e  —  07 

1 ,385e  -  02 

-7.738c  -  03 

1.342e  -  07 

1.043e  -  04 

- 1 .007e  -  04 

1.393e  -  05 

7.851e  -  06 

1.770e  -  02 

4.871c  -  04 

— 3.123e  -  05 

4.292c  -  04 

-2.549c-  03 

-3.550c  -  02 

2.258e  -  03 

-5.985c  -  05 

—  2.355e  -  03 

3.588c  -  02 

-2.273e  -  03 

—  5.48 1  -05 

1.184e  -  02 

—  2.322e  -  02 

Dl  = 

1 ,449e  -  03 

2.707e  -  04 

-1.925c  -  02 

-7.041e  -  03 

4.377e  -  04 

—  4 .396e  -  04 

-5.141c  -  03 

-7.797e  -  02 

4.806e  -  03 

-1.141c  -  04 

1 ,080e  -  02 

—  2 .006e  -  02 

1 .234c  -  03 

2.386e  -  04 

—  5.175e  -  04 

1.243c  -  01 

—  7  523e  -  03 

— 1.114e  -  05 

-1.916e  -  01 

-2.810c  -  05 

1.616e  -  06 

-3.771e  -  03 

-7.178c  -  05 

1.127e  -  01 

—  5.869e  -  03 

-1.181e  -  06 

—  1 .459e  -  03 

—  5.674e  -  03 

2  946e  -  04 

—  2.389e  -  05 

2.716e  -  01 

-2.096c  -  06 

1.052e  -  07 

4.242e  -  03 

— 9.120e  -  06 

— 4.698e  -  02 

2.175c  -  03 

- 1 .293c  -  07 

[  -1.770c -01 

2.299e  -  06 

-1.050e  -  07 

—  2.473c  -  03 

L  . 

j 

L. 

7.144e  -  07 

O.OOOe  +  00 

-3.946e  -  05 

9.562e 

-07 

O.OOOe  +  00 

-3.949c  -  05 

0.000c  +  00 

-2.170c  -08 

-4.200e  -  05 

O.OOOe  +  00 

-i.435e  -  08 

—  4 . 200e  -  05 

0.000c  +  00 

O.OOOe  +  00 

2.224e  -  01 

O.OOOe  +  00 

O.OOOe  +  00 

2.224e  -  01 

-1.564e -02 

2  570e  -  07 

—  2.994e  -  04 

—  1.710e 

-  02 

3.323e  -  07 

2  358e  -  0-1 

3  654e  -  07 

1.385e  -  02 

-3.112e  -  04 

3.967e 

-07 

1.778c  -  02 

2.829e  -  04 

1.043c  -  04 

-  1 ,007c  -  04 

-6.333c  -  02 

5.561e 

-  05 

— 6.235e  -  05 

7.231e  -  01 

1.770c  02 

4.871e  -  04 

5.976e  -  06 

-2.657e 

-  02 

—  4.646e  —  04 

—  7.02 le  -  04 

—  2.549e  -  03 

—  3.550e  -  02 

3.897c  -  04 

2.387e 

-  03 

3.479e  -  02 

—  4  811e  -  02 

—  2.355e  -  03 

3.588e  -  02 

1  220e  -  04 

2.399c 

-03 

-3.917c  -  02 

-1.555c  -  02 

1.184e  -  02 

—  2.322e  -  02 

3.80<:  04 

—  1  353e 

-  02 

1  884e  -  02 

—  5.238e  -  02 

—  1.925e  -  02 

—  7.041  e  -  03 

2.500e  -  05 

2.752e 

-  02 

5.873e  -  03 

—  3.534e  -  03 

—  5.141e  -  03 

—  7.797e  -  02 

—  2.456e  -  04 

4.962e 

-  03 

6.735c  -  02 

3.689c  -  02 

1  080e  -  02 

-2  006e  -  02 

3. 140e  -  04 

- 1 ,078c 

-  02 

1 ,835c  -  02 

-4.819c  -  02 

—  5.175e  -  04 

1  243e  -  01 

3.521e  -  05 

5.377e 

-  04 

—  8.600e  -  02 

-5.972c  -  03 

-1.916c  -  01 

-2.801c  -  05 

3.690e  -  05 

1  823e 

-  01 

7.747e  —  06 

-9.256e  -  03 

—  7.178e  -  05 

1.127e  -  01 

—  1.329e  -  04 

—  3.891  e 

-  05 

2  .829e  -  01 

9.886c  -  02 

-1  459e  -  03 

—  5.674e  -  03 

-2.518e  -  03 

-8.886c 

-  04 

—  1  464e  -  02 

’  937c  +  00 

2.716c  -  01 

-  2 ,096e  -  06 

-9.656e  -  06 

2.962e 

-  01 

-5.066e  -  06 

9.887c  -  03 

-9.120c  -  06 

—  4  .698e  -  02 

7.325e  -  06 

1.022c 

-  05 

3.637e  -  01 

—  6.133c  —  03 

-1.770e  -  01 


2.299e  —  06  -9.101c  -  06 


3.388e  -  01  -2.373e  -  05 


5.256e  -  03 
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1.417?  -  09 

3.503e  -  02 

3.468?  -  23 

1 .432?  -  09 

3.502?  - 

02 

1  663?  -  17 

3.807?  -  02 

O.OOOe  +  00 

-1.463?  -  24 

3.807?  -  02 

0.000?  +  00 

—  2.187e  -  18 

— 3.549e  -  07 

- 1 .698e  -  08 

-1.329?  -  26 

-3.546?  -  07 

-1.673?  - 

08 

1.440?  -  21 

—  1.522?  -  07 

—6.110?  -  03 

2.384?  -  15 

4.119?  -  06 

2.222?  - 

01 

3.055?  -  09 

—  T.738?  -  03 

1.342e  -  07 

-2.582?  -  15 

2.218?  -  01 

-5.184?  - 

06 

-3.310?  -  09 

1 .393?  -  05 

7.851?  -  06 

4.61  le  -  18 

-1.267?  -  03 

-1.195?  - 

03 

5.910?  -  12 

-3.123?  -  05 

4.292?  -  04 

1.387?  -  10 

7.993?  -  05 

1 .406?  - 

02 

1  778?  —  04 

2.25 8?  -  03 

—  5.985e  -  05 

-3.414?  -  11 

7.045?  -  03 

4.497?  - 

03 

-4.376?  -  05 

—  2.273e  -  03 

—  5.487e  -  05 

1.290?  -  10 

-3.180?  -  02 

3.039e  - 

03 

1.653?  -  04 

1.449?  -  03 

2.707e  -  04 

9.995?  -  12 

-1.190?  -  03 

-6.003?  - 

03 

1.274?  -05 

4.377e  -  04 

—  4 .396?  -  04 

O 

1 

o 

CN 

1.456?  -  03 

—  1.779e  - 

02 

2.588?  -04 

4.806?  -  03 

—  1.141e  -  04 

1.740?  -  10 

5.317?  -  02 

5.608?  - 

03 

2.230?  -  04 

1.234?  -  03 

2.386?  -  04 

8.846?  -  11 

2.096?  -  02 

-9.427?  - 

03 

1.134?  -  04 

—  7.523?  •  03 

-1.114?  -  05 

9.249?  -  11 

-7.424?  -  02 

1.910?  - 

04 

1.185?  -  04 

1.616?  -  06 

-3.771?  -  03 

-1.110?  -  11 

4.036?  -  05 

-1.311?  - 

01 

-1.423?  -  05 

—  5.869?  -  03 

—  1 . 18 1 e  -  06 

-7.828?  -  12 

1.680?  -  01 

-2.111?  - 

05 

- 1 .002?  -  05 

2.946?  -  04 

-2.389?  -  05 

4.332?  -  13 

-9.449?  -  03 

-3.073?  - 

04 

5.543?  —  07 

1.052?  -  07 

4.242?  -  03 

-1.050?  -  11 

-6.058?  -  06 

-1.686?  - 

01 

-1.343?  -  05 

2.175?  -  03 

-1.293?  -  07 

-5.344?  -  11 

1.906?  -  01 

-3.141?  - 

06 

-6.801?  -  05 

-1.050?  -  07 

-2.473?  -  03 

1 

00 

k) 

r# 

1 

-1.476?  -  05 

-1.763?  - 

01 

-1.051?  -  04 

The  cruciform  model  is  decentralized  into  x-axis  and  y-axis  systems  as  shown  in  Table 
D.l  [25], 


Table  D.l:  GHR  Decentralized  Model  Results. 


Model 

Inputs 

Outputs 

x-axis 

U i  .  U>2 

“2.  -5,  i'li  1J4,  VS-  !/ll 

y-axis 

U2,  Wi 

^1  •  ~4>  i/2'  i/5'  y 7'  yio 

APPENDIX  E 


CONTROLLER  PARAMETERS 


The  decentralized  second  order  controllers  which  are  applied  to  the  cruciform  structure 


are  given  in  the  following  form: 


’  H0 

G0  ' 

Ro 

O 

The  x-axis  controller  is  given  by 


—  5.5653c  —  01 

9.6470e 

-01 

— 8.2991e  -  01 

—  9.0911e  -  01 

3.2908e  - 

01 

—  1.1968e  +  01 

Cjr{s)  — 

— 4.3526e  -  01 

—  2.2147e 

-  01 

—  1.1 300e  +  00 

—  1 .2389c  +  00 

—  3.1525e  - 

01 

1.1474e  -t-  01 

-1.7181e  +  00 

-1.2874e  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 

The  y-axis  controller  is  given  by 


— 8.8423e 

-  01 

S.5382e  -  01 

— 3.60Sle  -  01 

-4.8280e  -  01 

9.2533e  -  01 

-2.6680e  +  01 

— 9.4732e 

-01 

—  1.3498e  -  01 

1.8310e  +  00 

2.537le  +  00 

9.2186e  -  02 

—  2.2745e  +  00 

— 9.7789e 

-02 

1.7635e  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 

O.OOOOe  +  00 
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